Longitudinal differential interferometric confocal microscopy

ABSTRACT

A differential interferometric confocal microscope for measuring an object, the microscope including a source-side pinhole array; a detector-side pinhole array; and an interferometer that images the array of pinholes of the source-side pinhole array onto a first array of spots located in front of an object plane located near where the object is positioned and onto a second array of spots behind the object plane, wherein the first and second arrays of spots are displaced from each other in both a direction normal to the object plane and a direction parallel to the object plane, the interferometer also imaging the first arrays of spots onto a first image plane that is behind the detector-side pinhole array and imaging the second array of spots onto a second image plane that is in front of the detector-side pinhole array wherein each spot of the imaged first array of spots is aligned with a corresponding different spot of the imaged second array of spots and a corresponding different pinhole of the detector-side pinhole array.

This application claims the benefit of U.S. Provisional Application No. 60/448,360, filed Feb. 19, 2003.

TECHNICAL FIELD

This invention relates to interferometric, confocal microscopes.

BACKGROUND OF THE INVENTION

There are a number of different forms of differential confocal microscopy. In one differential form, the Nomarski microscope measures one component of a conjugated quadratures of fields corresponding to the electrical interference signal of two images superimposed in an image plane. In another differential form, the conjugated quadratures of a dark field are measured one point at a time. In another differential form, the conjugated quadratures of each of two fields corresponding to two images are superimposed in an image plane one point at a time. In commonly owned U.S. Provisional Patent Application No. 60/447,254 (ZI-40) entitled “Transverse Differential Interferometric Confocal Microscopy” and U.S. patent application Ser. No. 10/778,371, filed Feb. 13, 2004 (ZI-40) also entitled “Transverse Differential Interferometric Confocal Microscopy” both of which are both by Henry A. Hill, it is taught how to practice transverse differential interferometric confocal microscopy. The contents of the cited U.S. Provisional Application and the U.S. Patent Application are herein incorporated in their entirety by reference.

However, neither the prior art nor the cited U.S. provisional patent application and cited U.S. Patent Application teach how to practice differential interferometric confocal microscopy wherein an array of conjugated quadratures of fields are measured jointly, where the components of each conjugated quadratures may be measured jointly, and where each conjugated quadratures represent a difference of conjugated quadratures of fields of converging beams subsequently scattered/reflected from a pair of locations on a substrate surface wherein one of the converging beams subsequently scattered/reflected or transmitted by the pair of locations is focused to an image plane located above the substrate surface and the second of the converging beams subsequently scattered/reflected from the pair of locations is focused to an image plane located below the substrate surface.

Also, prior art does not teach how to practice dark field differential interferometric confocal microscopy wherein an array of conjugated quadratures of fields are measured jointly, where the components of each conjugated quadratures may be measured jointly, where each conjugated quadratures represents a difference of conjugated quadratures of fields of converging beams subsequently scattered/reflected or transmitted by a pair of locations on a substrate surface wherein one of the converging beams subsequently scattered/reflected from the pair of locations is focused to an image plane located above the substrate surface and the second of the converging beams subsequently scattered/reflected or transmitted by the pair of locations is focused to an image plane located below the substrate surface, and where the nominal values of the conjugated quadratures of the array of conjugated quadratures is zero, i.e., the field that is being measured is nominally dark.

SUMMARY OF THE INVENTION

Embodiments of the present invention comprise interferometric confocal microscopy systems wherein an array of conjugated quadratures of fields are measured jointly, where the components of each conjugated quadratures may be measured jointly, where each conjugated quadratures represents a difference of conjugated quadratures of fields of converging beams subsequently scattered/reflected or transmitted by a pair of locations on a substrate surface wherein one of the converging beams subsequently scattered/reflected or transmitted by the pair of locations is focused to an image plane located above the substrate surface and the second of the converging beams subsequently scattered/reflected from the pair of locations is focused to an image plane located below the substrate surface, and where the nominal values of the conjugated quadratures of the array of conjugated quadratures may be adjusted as a set to be zero by controlling a single system parameter.

In general, in one aspect, the invention features a differential interferometric confocal microscope for measuring an object. The microscope includes a source-side pinhole array; a detector-side pinhole array; and an interferometer that images the array of pinholes of the source-side pinhole array onto a first array of spots located in front of an object plane located near where the object is positioned and onto a second array of spots behind the object plane, wherein the first and second arrays of spots are displaced from each other in both a direction normal to the object plane and a direction parallel to the object plane. The interferometer also images the first arrays of spots onto a first image plane that is behind the detector-side pinhole array and images the second array of spots onto a second image plane that is in front of the detector-side pinhole array wherein each spot of the imaged first array of spots is aligned with a corresponding different spot of the imaged second array of spots and a corresponding different pinhole of the detector-side pinhole array.

In general, in another aspect, the invention features a differential interferometric confocal microscope for measuring an object, wherein the microscope includes a source-side pinhole array; a detector-side pinhole array; and an interferometer that images each pinhole of the source-side pinhole array onto a corresponding different pair of two locations, one of which lies in a first object plane and the other of which lies in a second object plane that is parallel to and displaced from the first object plane, thereby generating a first image of the source-side pinhole array in the first object plane and a second image of the source-side pinhole array in the second object plane. The interferometer also projects a first array of return measurement beams from the first image and a second array of return measurement beams from the second image toward the detector-side pinhole array to produce a first array of converging beams and a second array of converging beams, wherein the detector-side pinhole array generates an array of conjugated quadratures of fields that is a difference of conjugated quadratures of fields of the first and second arrays of converging beams.

In general, in still another aspect, the invention features a differential interferometric confocal microscope for measuring an object and which has, in the vicinity of where the object being measured is to be located, a first object plane and a second object plane that is displaced from and parallel to the first object plane. The microscope includes: a source-side pinhole array; a detector-side pinhole array; and an interferometer that receives a beam from a selected pinhole of the source-side pinhole array and converges a first part of that received beam onto a corresponding first location in the first object plane and a second part of that received beam onto a corresponding second location in the second object plane. The interferometer is further arranged to receive a first return beam from the first location and a second return beam from the second location and converge at least a part of each of the first and second return beams onto a corresponding pinhole of the detector-side pinhole array to produce a difference of conjugated quadratures of fields of the first and second return beams converging on that corresponding pinhole, wherein the selected pinhole is any pinhole of the source-side pinhole array.

In general, in still yet another aspect, the invention features a differential interferometric confocal microscope for measuring an object. The microscope includes: a source-side pinhole array for producing an array of input beams; a detector-side pinhole array; and an interferometer. The interferometer includes a first optical element providing a first reflecting surface; a second optical element providing a second reflecting surface; and a beam splitter positioned between the first and second optical elements, wherein the beam splitter produces from the array of input beams a first array of measurement beams and a second array of measurement beams, wherein the first reflecting surface participates in focusing the first array of measurement beams onto a first array of locations on a first object plane in object space and the second reflecting surface participates in focusing the second array of measurement beams onto a second array of locations on a second object plane in object space, the first and second object planes being parallel to and displaced from each other. Also, the first array of measurement beams generates a first array of return beams from the object and the second array of measurement beams generates a second array of return beams from the object, and the first reflecting element participates in producing from the first array of return beams a first array of converging beams that converge to a first array of spots on a first image plane and the second reflecting element participates in producing from the second array of return beams a second array of converging beams that converge onto a second array of spots on a second image plane. The first and second image planes are adjacent to and on opposite sides of the detector-side pinhole array, and wherein the detector-side pinhole array combines the first and second arrays of converging beams to form an array of output beams.

Other embodiments include one or more of the following features. A single pinhole array serves as both the source-side pinhole array and the detector-side pinhole array. The first optical element is located between the single pinhole array and the beam splitter and wherein the second optical element is located between a location at which the object is positioned during use and the beam splitter. The first reflecting surface has a center of curvature for which there is a corresponding conjugate as viewed through the beam splitter, and the second reflecting surface has a center of curvature that is displaced relative to the corresponding conjugate of the center of curvature of the first reflecting surface. More specifically, the conjugate of the center of curvature of the first reflecting surface and the center of curvature of the second reflecting surface are displaced from each other in a first direction that is normal to a plane defined by the beam splitter and in a second direction that is parallel to the plane defined by the beam splitter. The reflecting surface participates in focusing the first array of measurement beams via the beam splitter onto the first array of locations and the second reflecting surface participates in focusing the second array of measurement beams via the beam splitter onto the second array of locations. The first reflecting element participates in combination with the beam splitter in producing the first array of converging beams and the second reflecting element participates in combination with the beam splitter in producing the second array of converging beams. The first reflecting surface is substantially concentric with a point on the object. The second optical element provides a refracting surface positioned between the object and the beam splitter to receive light rays from the object. The first reflecting surface substantially conforms to a sphere having a first radius and the refracting surface conforms to a sphere having a second radius, wherein the first radius is greater than the second radius. The first optical element provides a refracting surface positioned between the beam splitter and the single pinhole array. The second reflecting surface is substantially concentric with an image point on the single pinhole array. The second reflecting surface substantially conforms to a sphere having a first radius and the refracting surface conforms to a sphere having a second radius, wherein the first radius is greater than the second radius. The single pinhole array is a two-dimensional array; more specifically, a two-dimensional array of equally-spaced holes, which are circular apertures. The first and second object planes are separated from each other on the order of the longitudinal resolution of the differential confocal interferometric microscope.

An advantage of at least one embodiment of the present invention is that the fields of beams scattered/reflected or transmitted by a pair of locations on a substrate surface are generated by a single confocal pinhole.

Another advantage of at least one embodiment of the present invention is that reference beam components of an array of reference beams used in generation of electrical interference signals corresponding to measured conjugated quadratures of fields of beams scattered/reflected or transmitted by a pair of locations on a substrate are identical.

Another advantage of at least one embodiment of the present invention is that components of background beams generated by measurement beam components subsequently scattered/reflected or transmitted at a pair of locations on a substrate surface are substantially identical at a confocal pinhole.

Another advantage of at least one embodiment of the present invention is that the spatial filtering of fields of beams scattered/reflected or transmitted at a pair of locations on a substrate surface is performed by a single confocal pinhole.

Another advantage of at least one embodiment of the present invention is that information about a substrate surface is obtained with an interferometric confocal imaging system operating in a dark field mode.

Another advantage of at least one embodiment of the present invention is that information about a substrate surface is obtained with reduced systematic and statistical errors.

Another advantage of at least one embodiment of the present invention is the generation of a significant increase in throughput because the intensity of an input beam may be significantly increased without saturation of a detector system.

Another advantage of at least one embodiment of the present invention is that an array of conjugated quadratures of the fields of beams scattered/reflected or transmitted by a pair of locations on a substrate surface is measured jointly and the components of each conjugated quadratures may be measured jointly.

Another advantage of at least one embodiment of the present invention is that information is obtained about critical dimensions and locations of sub-wavelength artifacts on a substrate surface.

Another advantage of at least one embodiment of the present invention is that information is obtained about the sizes and locations of sub-wavelength defects on a substrate surface.

Another advantage of at least one embodiment of the present invention is that information may be obtained about a longitudinal derivative of a profile of a substrate surface.

Another advantage of at least one embodiment of the present invention is that information may be obtained about one-dimensional and two-dimensional profiles of a substrate surface.

Another advantage of at least one embodiment of the present invention is that imaging of a longitudinal gradient of a substrate surface profile with a lateral resolution of the order of 100 nm and a longitudinal resolution of the order of 200 nm may be obtained with a working distance of the order of a mm.

The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 a is a diagram of an interferometric system used to make differential measurements of conjugated quadratures of fields of beams scattered/reflected or transmitted by a substrate.

FIG. 1 b is a schematic diagram of a beam-conditioner configured to operate in a two-frequency generator and phase-shifter.

FIG. 1 c is a schematic diagram of a beam-conditioner configured to operate in a two-frequency generator and frequency-shifter.

FIG. 2 a is a schematic diagram of a confocal microscope system.

FIG. 2 b is a schematic diagram of catadioptric imaging system.

FIG. 2 c is a schematic diagram of beams focused to spots at a pinhole array used in a confocal microscope system.

FIG. 2 d is a schematic diagram of beams focused to spots in a catadioptric imaging system.

FIG. 2 e is a schematic diagram of beams focused to spots in a catadioptric imaging system.

FIG. 2 f is a schematic diagram of beams focused to spots in an image plane.

FIG. 3 is a schematic diagram of an interferometric confocal imaging system used to make differential measurements of conjugated quadratures of fields of beams scattered/reflected by a substrate.

DETAILED DESCRIPTION

In general, an array of conjugated quadratures of fields is measured interferometrically by a differential confocal interferometer and detector system wherein each conjugated quadratures comprises a difference of conjugated quadratures of fields of converging beams scattered/reflected or transmitted by a pair of locations on a substrate surface. One of the converging beams subsequently scattered/reflected from the pair of locations is converging to a spot in a first image plane located above the substrate surface and the other of the converging beams subsequently scattered/reflected from the pair of locations is converging to a spot in a second image plane located below the substrate surface. The array of conjugated quadratures is measured jointly, i.e., simultaneously, and the components of each conjugated quadratures may be measured jointly. The separation of the first and second image planes is of the order of the longitudinal resolution of the differential confocal interferometer and detector system. Each pair of locations generally has a relative displacement of the order of the transverse resolution of the differential confocal interferometer and detector system in a direction nominally tangent to the substrate surface. The relative phases of the converging beams subsequently scattered/reflected by the substrate surface may be adjusted as a set by control of a single system parameter so that the conjugated quadratures of the array of conjugated quadratures are nominally zero, i.e., information may be obtained about the substrate surface with the interferometer and detector system operating in a dark field mode. Operation in a dark field mode leads to both reduced systematic and statistical errors in the information and to increased throughput. The information may comprise one-dimensional and two-dimensional longitudinal profiles of a substrate surface, one-dimensional and two-dimensional longitudinal derivatives of a profile of a substrate surface; critical dimensions of features or artifacts on a substrate surface; and the size and location of sub-wavelength defects on a substrate surface.

In one embodiment, an image plane of an interferometric confocal imaging system comprises a superposition of two images of a substrate surface wherein each of the two superimposed images corresponds to a partially defocused image of a location on the substrate surface and the two partially defocused images are displaced transversely relative to each other. An array of conjugated quadratures of fields representing the superimposed images are measured jointly and the components of each conjugated quadratures may be measured jointly. Each pair of locations on the substrate surface corresponding to a conjugated quadratures of the array of conjugated quadratures generally has a relative transverse displacement of the order of 3 times the size of the transverse resolution of the interferometric confocal imaging system in a direction nominally tangent to the surface of the imaged section. The longitudinal separation of the image planes of the beams generating the partially defocused images of the substrate surface is of the order of the longitudinal resolution of the interferometric confocal imaging system. The respective conjugated quadrature of a field is a sin φ when the quadrature x(φ) of the field is a cos φ.

In another embodiment, the relative phases of beams subsequently scattered/reflected or transmitted at a pair of locations on a substrate surface may be adjusted by a single confocal imaging system parameter so that conjugated quadratures of an array of conjugated quadratures of fields of the scattered/reflected or transmitted beams from the respective pair of locations are nominally zero, i.e., information is obtained about the substrate surface with the interferometric imaging system operating in a dark field mode.

Operation in a dark field mode leads to both reduced systematic and statistical errors in the information. When operating in a dark field mode, a measured conjugated quadratures of fields of beams scattered/reflected or transmitted from a respective pair of locations comprising a sub-wavelength artifact in an otherwise locally isotropic substrate surface represents information only about the sub-wavelength artifact relative to a reference sub-wavelength artifact. The reference artifact has properties of the otherwise locally isotropic substrate surface and dimensions similar to those of the artifact. Accordingly, properties measured include information about critical dimensions and location of the sub-wavelength artifact on the substrate.

Also when operating in the dark field mode, a measured conjugated quadratures of fields of beams scattered/reflected or transmitted by a respective pair of locations on a substrate surface comprising a sub-wavelength defect in an otherwise locally isotropic substrate surface represents information only about the sub-wavelength defect relative to a reference sub-wavelength defect. The reference defect has properties of the otherwise locally isotropic substrate surface and dimensions similar to those of the defect. Accordingly, properties measured include information about dimensions and location of the sub-wavelength defect on the substrate.

A general description of embodiments incorporating aspects of the present invention will first be given wherein the embodiments comprise an interferometer system that uses either a single-, double, bi-, or quad-homodyne detection method and a partially defocused image of an array of a first a set of locations on a substrate surface and a partially defocused image of an array of a second set of locations on the substrate surface are superimposed on an image plane of the interferometer system. The transverse spacing between pairs of locations comprising corresponding locations of the first and second sets of locations on the substrate surface is of the order of at least 3 times the size of the transverse resolution of the interferometer system. The longitudinal separation of the image planes of the beams generating the partially defocused images of the substrate surface is of the order of the longitudinal resolution of the interferometer system. There is a one-to-two mapping of a location in the superimposed image space to the two locations of pair of locations.

Referring to FIG. 1 a, an interferometer system is shown diagrammatically comprising an interferometer generally shown as numeral 10, a source 18, a beam-conditioner 22, detector 70, an electronic processor and controller 80, and a measurement object or substrate 60. Source 18 and beam conditioner 22 generate input beam 24 comprising one or more frequency components. Source 18 is a pulsed source. Two or more of the frequency components of input beam 24 may be coextensive in space and may have the same temporal window function.

Reference and measurement beams are generated in interferometer 10 for each of the frequency components of beam 24. The measurement beam generated in interferometer 10 is one component of beam 28 and imaged to form a partially defocused image on the surface of substrate 60 to form an array of pairs of partially defocused images. Beam 28 further comprises a return reflected/scattered measurement beam that is generated by the reflection/scattering or transmission of the measurement beam component of beam 28 at the array of pairs of partially defocused images on the surface of substrate 60. Interferometer 10 superimposes the two arrays of components of the return measurement beam corresponding to the two arrays components of beam 28 reflected/scattered or transmitted at the arrays of the pairs of defocused images to form a single array of superimposed images of return measurement beam components of beam 28. The return measurement beam components of beam 28 are subsequently combined with the reference beam in interferometer 10 to form output beam 32.

Output beam 32 is detected by detector 70 to generate an electrical interference signal 72. Detector 70 may comprise an analyzer to select common polarization states of the reference and return measurement beam components of beam 32 to form a mixed beam. Alternatively, interferometer 10 may comprise an analyzer to select common polarization states of the reference and return measurement beam components such that beam 32 is a mixed beam.

Two different modes are described for the acquisition of the electrical interference signals 72. The first mode to be described is a step and stare mode wherein substrate 60 is stepped between fixed locations corresponding to locations where image information is desired. The second mode is a scanning mode. In the step and stare mode for generating a one-dimensional and a two-dimensional surface profile of substrate 60, substrate 60 mounted in wafer chuck 84/stage 90 is translated by stage 90. The position of stage 90 is controlled by transducer 82 according to servo control signal 78 from electronic processor and controller 80. The position of stage 90 is measured by metrology system 88 and position information acquired by metrology system 88 is transmitted to electronic processor and controller 80 to generate an error signal for use in the position control of stage 90. Metrology system 88 may comprise for example linear displacement and angular displacement interferometers and cap gauges.

Electronic processor and controller 80 translates wafer stage 90 to a desired position and then acquires a set of four electrical interference signal values. After the acquisition of the sequence of four electrical interference signals, electronic processor and controller 80 then repeats the procedure for the next desired position of stage 90. The elevation and angular orientation of substrate 60 is controlled by transducers 86A and 86B.

The second mode for the acquisition of the electrical interference signal values is next described wherein the electrical interference signal values are obtained with the position of stage 90 scanned in one or more directions. In the scanning mode, source 18 is pulsed at times controlled by signal 92 from signal processor and controller 80. Source 18 is pulsed at times corresponding to the registration of the conjugate image of confocal pinholes or pixels of detector 70 with positions on and/or in substrate 60 for which image information is desired.

There will be a restriction on the duration or “pulse width” of a beam pulse τ_(p1) produced by source 18 as a result of the continuous scanning mode used in the third variant of the first embodiment. Pulse width τ_(p1) will be a parameter that in part controls the limiting value for spatial resolution in the direction of a scan to a lower bound of τ_(p1)V,  (1) where V is the scan speed. For example, with a value of τ_(p1)=50 nsec and a scan speed of V=0.20 m/sec, the limiting value of the spatial resolution τ_(p1)V in the direction of scan will be τ_(p1)V=10 nm.  (2)

Pulse width τ_(p1) will also determine the minimum frequency difference that can be used in the bi- and quad-homodyne detection methods. In order that there be no contributions to the electrical interference signals from interference between fields of conjugated quadratures, the minimum frequency spacing Δf_(min) is expressed as

$\begin{matrix} {{\Delta\; f_{\min}}\operatorname{>>}{\frac{1}{\tau_{p1}}.}} & (3) \end{matrix}$ For an example of τ_(p1)=50 nsec, 1/τ_(p1)=20 MHz.

For certain embodiments, the frequencies of input beam 24 are controlled by signals 74 and/or 92 from signal processor and controller 80 to correspond to the frequencies that will yield the desired phase shifts between the reference and return measurement beam components of output beam 32. Alternatively in certain other embodiments, the relative phases of reference and measurement beam components of input beam 24 are controlled by signal 74 and/or 92 from signal processor and controller 80 to correspond to the desired phase shifts between the reference and return measurement beam components of output beam 32. In the first mode, i.e., the step and stare mode, each set of the sets of arrays of four electrical interference signal values corresponding to the set of four phase shift values are generated by a single pixel of detector 70 for single- and bi-homodyne detection method, by two pixels of detector 70 for the quad-homodyne detection method, and by four pixels of detector 70 for the double-homodyne detection methods. In the second mode for the acquisition of the electrical interference signal values, each corresponding set of four electrical interference signal values are generated by a conjugate set of four different pixels of detector 70 for each of the four homodyne detection methods. Thus in the first mode of acquisition, the differences in pixel efficiency are compensated in the signal processing by signal processor and controller 80 for the double-, bi-, and quad-homodyne detection methods and in the second mode of acquisition, the differences in pixel efficiency and the differences in sizes of pinholes in confocal pinhole arrays are compensated in the signal processing by signal processor and controller 80 as described in the subsequent description of the homodyne detection methods. The joint measurements of conjugated quadratures of fields generated by electric processor and controller 80 are subsequently described in the description of the bi- and quad-homodyne detection methods.

In practice, known phase shifts are introduced between the reference and measurement beam components of output beam 32 by two different techniques. In one technique, phase shifts are introduced between the reference and measurement beam components for each of the at least two frequency components by source 18 and beam-conditioner 22 as controlled by signals 92 and 74, respectively, from electronic processor and controller 80. In the second technique, phase shifts are introduced between the reference and measurement beam components for each of the at least two frequency components as a consequence of frequency shifts introduced to the frequency components of input beam 24 by source 18 and beam-conditioner 22 as controlled by signals 92 and 74, respectively, from electronic processor and controller 80.

There are different ways to configure source 18 and beam-conditioner 22 to meet the input beam requirements of different embodiments. Reference is made to FIG. 1 b where beam-conditioner 22 is configured as a two-frequency generator and a phase-shifter and source 18 is configured to generate beam 20 with one frequency component. The two-frequency generator and phase-shifter configuration comprises acousto-optic modulators 1020, 1026, 1064, and 1068; polarizing beam-splitters 1030, 1042, 1044, and 1056; phase-shifters 1040 and 1052; half wave phase retardation plates 1072 and 1074; non-polarizing beam-splitter 1070, and mirrors 1036, 1038, 1050, 1054, and 1056.

Input beam 20 is incident on acousto-optic modulator 1020 with a plane of polarization parallel to the plane of FIG. 1 b. A first portion of beam 20 is diffracted by acousto-optic modulator 1020 as beam 1022 and then by acousto-optic modulator 1026 as beam 1028 having a polarization parallel to the plane of FIG. 1 b. A second portion of beam 20 is transmitted as a non-diffracted beam 1024 having a plane of polarization parallel to the plane of FIG. 1 b. The acoustic power to acousto-optic modulator 1020 is adjusted such that beams 1022 and 1024 have nominally the same intensity.

Acousto-optic modulators 1020 and 1026 may be of either the non-isotropic Bragg diffraction type or of the isotropic Bragg diffraction type. The frequency shifts introduced by acousto-optic modulators 1020 and 1026 are of the same sign and equal to ¼ of the desired frequency shift between the two frequency components of input beam 24. Also the direction of propagation of beam 1028 is parallel to the direction of propagation of beam 1024.

Beam 1024 is diffracted by acousto-optic modulators 1064 and 1068 as beam 1082 having a polarization parallel to the plane of FIG. 1 b. Acousto-optic modulators 1064 and 1068 may be of either the non-isotropic Bragg diffraction type or of the isotropic Bragg diffraction type. The frequency shifts introduced by acousto-optic modulators 1064 and 1068 are of the same sign and equal to ¼ of the desired frequency shift between the two frequency components of input beam 24. Also the direction of propagation of beam 1082 is parallel to the direction of propagation of beam 1024.

Beams 1028 and 1082 are incident on half-wave phase retardation plates 1072 and 1074, respectively, and transmitted as beams 1076 and 1078, respectively. Half-wave phase retardation plates 1072 and 1074 are oriented such that the planes of polarization of beams 1076 and 1078 are at 45 degrees to the plane of FIG. 1 b. The components of beams 1076 and 1078 polarized parallel to the plane of FIG. 1 b will be used as the measurement beam components in interferometer 10 and the components of beams 1076 and 1078 polarized orthogonal to the plane of FIG. 1 b will be used as the reference beam components in interferometer 10.

Continuing with reference to FIG. 1 b, beam 1076 is incident on polarizing beam-splitter 1044 and the respective measurement and reference beam components transmitted and reflected, respectively, as beams 1046 and 1048, respectively. Measurement beam component 1046 is transmitted by polarizing beam-splitter 1056 as a measurement beam component of beam 1058 after reflection by mirror 1054. Reference beam component 1048 is reflected by polarizing beam-splitter 1056 as reference beam component of beam 1058 after reflection by mirror 1050 and transmission by phase-shifter 1052. Beam 1058 is incident on beam-splitter 1070 and a portion thereof is reflected as a component of beam 24.

Beam 1078 is incident on polarizing beam-splitter 1030 and the respective measurement and reference beam components transmitted and reflected, respectively, as beams 1032 and 1034, respectively. Measurement beam component 1032 is transmitted by polarizing beam-splitter 1042 as a measurement beam component of beam 1060 after reflection by mirror 1036. Reference beam component 1034 is reflected by polarizing beam splitter 1042 as reference beam component of beam 1060 after reflection by mirror 1038 and transmission by phase-shifter 1040. Beam 1060 is incident on beam-splitter 1070 and a portion thereof is transmitted as a component of beam 24 after reflection by mirror 1056.

Phase-shifters 1052 and 1040 introduce phase shifts between respective reference and measurement beams according to signal 74 from electronic processor and controller 80 (see FIG. 1 a). The schedule of the respective phase shifts is described in the subsequent discussions of homodyne detection methods. Phase-shifters 1052 and 1040 may be for example of the optical-mechanical type comprising for example prisms and piezoelectric translators or of the electro-optical modulator type.

Beam 24 that exits the two-frequency generator and phase shift configuration of beam-conditioner 22 comprises one reference beam and measurement beam having one frequency, a second reference beam and measurement beam having a second frequency component, and relative phases of the reference beams and the measurement beams that are controlled by electronic processor and controller 80.

Continuing with a description of different ways to configure source 18 and beam-conditioner 22 to meet the input beam requirements of different embodiments, reference is made to FIG. 1 c where beam-conditioner 22 is configured as a two-frequency generator and a frequency shifter. The two-frequency generator and frequency-shifter configuration comprises acousto-optic modulators 1120, 1126, 1130, 1132, 1142, 1146, 1150, 1154, 1058, and 1062; beam-splitter 1168; and mirror 1166.

Source 18 is configured to generate beam 20 with a single frequency component. Beam 20 is incident on acousto-optic modulator 1120 with a plane of polarization parallel to the plane of FIG. 1 c. A first portion of beam 20 is diffracted by acousto-optic modulator 1120 as beam 1122 and then by acousto-optic modulator 1126 as beam 1128 having a polarization parallel to the plane of FIG. 1 c. A second portion of beam 20 is transmitted as a non-diffracted beam 1124 having a plane of polarization parallel to the plane of FIG. 1 c. The acoustic power to acousto-optic modulator 1120 is adjusted such that beams 1122 and 1124 have nominally the same intensity.

Acousto-optic modulators 1120 and 1126 may be of either the non-isotropic Bragg diffraction type or of the isotropic Bragg diffraction type. The frequency shifts introduced by acousto-optic modulators 1120 and 1126 are of the same sign and equal to ½ of a frequency shift Δf that will generate in interferometer 10 a relative π/2 phase shift between a corresponding reference beam and a measurement beam that have a relative change in frequency equal to the frequency shift. The direction of propagation of beam 1128 is parallel to the direction of propagation of beam 1124.

Continuing with FIG. 1 c, beam 1128 is incident on acousto-optic modulator 1132 and is either diffracted by acousto-optic modulator 1132 as beam 1134 or transmitted by acousto-optic modulator 1132 as beam 1136 according to control signal 74 (see FIG. 1 a) from electronic processor and controller 80. When beam 1134 is generated, beam 1134 is diffracted by acousto-optic modulators 1142, 1146, and 1150 as a frequency-shifted beam component of beam 1152. The frequency shifts introduced by acousto-optic modulators 1132, 1142, 1146, and 1150 are all in the same direction and equal in magnitude to Δf/2. Thus the net frequency shift introduced by acousto-optic modulators 1132, 1142, 1146, and 1150 is ±2Δf and will generate a relative π phase between the respective reference and measurement beams in interferometer 10. The net frequency shift introduced by acousto-optic modulators 1120, 1126, 1132, 1142, 1146, and 1150 is Δf±2Δf and will generate a respective relative phase shift of π/2±π between the respective reference and measurement beams in interferometer 10.

When beam 1136 is generated, beam 1136 is transmitted by acousto-optic modulator 1150 according to control signal 74 from electronic processor and controller 80 as a non-frequency shifted beam component of beam 1152 with respect to beam 1128. The frequency shift introduced by acousto-optic modulators 1120, 1126, and 1150 is Δf and will generate a respective relative phase shift of π/2 between the respective reference and measurement beams in interferometer 10.

Beam 1124 is incident on acousto-optic modulator 1130 and is either diffracted by acousto-optic modulator 1130 as beam 1140 or transmitted by acousto-optic modulator 1130 as beam 1138 according to control signal 74 from electronic processor and controller 80. When beam 1140 is generated, beam 1140 is diffracted by acousto-optic modulators 1154, 1158, and 1162 as a frequency-shifted beam component of beam 1164. The frequency shifts introduced by acousto-optic modulators 1130, 1154, 1158, and 1162 are all in the same direction and equal to ±Δf/2. Thus the net frequency shift introduced by acousto-optic modulators 1130, 1154, 1158, and 1162 is ±Δf/2 and will generate a relative phase shift of π between the respective reference and measurement beams on transit through interferometer 10. The net frequency shift introduced by acousto-optic modulators 1120, 1130, 1154, 1158, and 1162 is ±Δf/2 and will generate a respective relative phase shift of ±π between the respective reference and measurement beams on transit through interferometer 10

When beam 1138 is generated, beam 1138 is transmitted by acousto-optic modulator 1162 according to control signal 74 from electronic processor and controller 80 as a non-frequency shifted beam component of beam 1164. The frequency shift introduced by acousto-optic modulators 1120, 1130, and 1162 is 0 and will generate a respective relative phase shift of 0 between the respective reference and measurement beams on transit through interferometer 10.

Beams 1152 and 1164 may be used directly as input beam 24 when an embodiment requires spatially separated reference and measurement beams for an input beam. When an embodiment requires coextensive reference and measurement beams as an input beam, beam 1152 and 1164 are next combined by beam-splitter 1168 to form beam 24. Acousto-optic modulators 1120, 1126, 1130, 1132, 1142, 1146, 1150, 1154, 1058, and 1062 may be either of the non-isotropic Bragg diffraction type or of the isotropic Bragg diffraction type. Beams 1152 and 1164 are both polarized in the plane of FIG. 1 c for either non-isotropic Bragg diffraction type or of the isotropic Bragg diffraction type and beam-splitter 1168 is of the non-polarizing type.

With a continuation of the description of different ways to configure source 18 and beam-conditioner 22 to meet the input beam requirements of different embodiments, source 18 will preferably comprise a pulsed source. There are a number of different ways for producing a pulsed source [see Chapter 11 entitled “Lasers”, Handbook of Optics, 1, 1995 (McGraw-Hill, New York) by W. Silfvast]. Each pulse of source 18 may comprise a single pulse or a train of pulses such as generated by a mode locked Q-switched Nd:YAG laser. A single pulse train is referenced herein as a pulse sequence and a pulse and a pulse sequence are used herein interchangeably.

Source 18 may be configured in certain embodiments to generate one or more frequencies by techniques such as described in a review article entitled “Tunable, Coherent Sources For High-Resolution VUV and XUV Spectroscopy” by B. P. Stoicheff, J. R. Banic, P. Herman, W. Jamroz, P. E. LaRocque, and R. H. Lipson in Laser Techniques for Extreme Ultraviolet Spectroscopy, T. J. McIlrath and R. R. Freeman, Eds., (American Institute of Physics) p 19 (1982) and references therein. The techniques include for example second and third harmonic generation and parametric generation such as described in the articles entitled “Generation of Ultraviolet and Vacuum Ultraviolet Radiation” by S. E. Harris, J. F. Young, A. H. Kung, D. M. Bloom, and G. C. Bjorklund in Laser Spectroscopy I, R. G. Brewer and A. Mooradi, Eds. (Plenum Press, New York) p 59, (1974) and “Generation of Tunable Picosecond VUV Radiation” by A. H. Kung, Appl. Phys. Lett. 25, p 653 (1974). The contents of the three cited articles are herein incorporated in their entirety by reference.

The output beams from source 18 comprising two or four frequency components may be combined in beam-conditioner 22 by beam-splitters to form coextensive measurement and reference beams that are either spatially separated or coextensive as required in various embodiments. When source 18 is configured to furnish two or four frequency components, the frequency shifting of the various components required in certain embodiments may be introduced in source 18 for example by frequency modulation of input beams to parametric generators and the phase shifting of reference beams relative to measurement beams in beam-conditioner 22 may be achieved by phase shifters of the optical-mechanical type comprising for example prisms or mirrors and piezoelectric translators or of the electro-optical modulator type.

The general description of embodiments incorporating various aspects of the present invention is continued with reference to FIG. 1 a. Input beam 24 is incident on interferometer 10 wherein reference beams and measurement beams are present in input beam 24 or are generated from input beam 24 in interferometer 10. The reference beams and measurement beams comprise two arrays of reference beams and two arrays of measurement beams wherein the arrays may comprise arrays of one element. The arrays of measurement beams are incident on or focused on and/or in substrate 60 and arrays of return measurement beams are generated by reflection/scattering and/or transmission by the substrate. In the case of single element arrays for the reference beams and measurement beams, the measurement beams are generally reflected or transmitted by substrate 60. The arrays of reference beams and return measurement beams are combined by a beam-splitter to form two arrays of output beam components. The arrays of output beam components are mixed with respect to state of polarization either in interferometer 10 or in detector 70. The arrays of output beams are subsequently focused to spots on pixels of a multi-pixel or single pixel detector as required and detected to generate electrical interference signal 72.

There are four different implementations of the homodyne detection method that are used in interferometric embodiments. The four different implementations are referred to as single-, double-, bi-, and quad-homodyne detection methods. For the single-homodyne detection method, input beam 24 comprises a single frequency component and a set of four measurements of the array of electrical interference signals 72 is made. For each of the four measurements of the array of electrical interference signals 72, a known phase shift is introduced between the reference beam component and respective return measurement beam components of output beam 32. The subsequent data processing procedure used to extract the conjugated quadratures of the reflected and/or scattered or transmitted return measurement beam for an input beam comprising a single frequency component is described for example in commonly owned U.S. Pat. No. 6,445,453 (ZI-14) entitled “Scanning Interferometric Near-Field Confocal Microscopy” by Henry A. Hill, the contents of which are herein incorporated in their entirety by reference.

The double-homodyne detection method uses input beam 24 comprising four frequency components and four detectors to obtain measurements of electrical interference signals that are subsequently used to obtain conjugated quadratures. Each detector element of the four detector elements obtains a different one of the four electrical interference signal values with the four electrical interference signal values obtained simultaneously to compute the conjugated quadratures for a field. Each of the four electrical interference signal values contains only information relevant to one orthogonal component of the conjugated quadratures. The double-homodyne detection used herein is related to the detection methods such as described in Section IV of the article by G. M D'ariano and M G. A. Paris entitled. “Lower Bounds On Phase Sensitivity In Ideal And Feasible Measurements,” Phys. Rev. A 49, 3022–3036 (1994). Accordingly, the double-homodyne detection method does not make joint determinations of conjugated quadratures of fields wherein each electrical interference signal value contains information simultaneously about each of two orthogonal components of the conjugated quadratures.

The bi- and quad-homodyne detection methods obtain measurements of electrical interference signals wherein each measured value of an electrical interference signal contains simultaneously information about two orthogonal components of conjugated quadratures. The two orthogonal components correspond to orthogonal components of conjugated quadratures such as described in cited U.S Provisional Patent Application No. 60/442,858 (ZI-47) and cited U.S. patent application Ser. No. 10/765,369, filed Jan. 27, 2004 (ZI-47) entitled “Apparatus and Method for Joint Measurements of Conjugated Quadratures of Fields of Reflected/Scattered and Transmitted Beams by an Object in Interferometry.”

Conjugated quadratures of fields of the return measurement beam are obtained by single-, double-, bi-, and quad- homodyne detection methods in the interferometric embodiments. For each of the homodyne detection methods, a set of four measurements of the array of electrical interference signals 72 is made. For each of the four measurements of the array of electrical interference signals 72, a known phase shift is introduced between the reference beam components and respective return measurement beam components of output beam 32. A nonlimiting example of a known set of phase shifts comprise 0, π/4, π/2, and 3π/2 radians, mod 2π.

Input beam 24 comprises for interferometric embodiments one frequency component for the single-homodyne detection method. For the bi-homodyne detection method, input beam 24 comprises two frequency components and for double- and quad-homodyne detection methods, input beam 24 comprises four frequency components. The phase shifts are generated by either shifting the frequencies of frequency components of input beam 24 between known frequency values or by introducing phase shifts between the reference and measurement beam components of input beam 24. In certain of the interferometric embodiments, there is a difference between the optical path lengths of the reference beam components and the respective return beam components of output beam components such for output beam 32 in interferometer 10. As a consequence, a change in frequency of a frequency component of input beam 24 will generate a relative phase shift between the corresponding reference beam components and the respective return beam components of output beam 32.

For an optical path difference L between the reference beam components and the respective return measurement beam components of output beam 32, there will be for a frequency shift Δf a corresponding phase shift Φ where

$\begin{matrix} {\varphi = {2\;\pi\;{L\left( \frac{\Delta\; f}{c} \right)}}} & (4) \end{matrix}$ and c is the free space speed of light. Note that L is not a physical path length difference and depends for example on the average index of refraction of the measurement beam and the return measurement beam paths. For an example of a phase shift Φ=π, 3π, 5π, . . . and a value of L=0.25 m, the corresponding frequency shifts are Δf=600 MHz, 1.8 GHz, 3.0 GHz, . . . .

The frequencies of components of input beam 24 are determined by the mode of operation of source 18 and of beam-conditioner 22 according to control signals 92 and 74, respectively, generated by electronic processor and controller 80.

Referring to the bi-homodyne detection method used in some embodiments, a set of four electrical interference signal values are obtained for each pair of spots in or on substrate 60 being imaged such as described in commonly owned U.S. Provisional Patent Application No. 60/442,858 (ZI-47) and entitled “Apparatus and Method for Joint Measurements of Conjugated Quadratures of Fields of Reflected/Scattered Beams by an Object in Interferometry” and U.S. patent application Ser. No. 10/765,369, filed Jan. 27, 2004 (ZI-47) and entitled “Apparatus and Method for Joint Measurements of Conjugated Quadratures of Fields of Reflected/Scattered and Transmitted Beams by an Object in Interferometry” both of which are by Henry A. Hill. The contents of both the cited U.S. Provisional Patent Application and the U.S. Patent Application are herein incorporated in their entirety by reference. The set of four electrical interference signal values S_(j), j=1,2,3,4 used for obtaining conjugated quadratures of fields for a single a spot on and/or in a substrate being imaged is represented for the bi-homodyne detection within a scale factor by the formula

$\begin{matrix} {S_{j} = {P_{j}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{j}^{2}{A_{1}}^{2}} + {\zeta_{j}^{2}{B_{1}}^{2}} + {\eta_{j}^{2}{C_{1}}^{2}} + {\zeta_{j}\eta_{j}2{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}} +} \\ {{\xi_{j}\zeta_{j}2{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}} + {ɛ_{j}\xi_{j}\eta_{j}2{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \end{matrix} \\ {{\xi_{j}^{2}{A_{2}}^{2}} + {\zeta_{j}^{2}{B_{2}}^{2}} + {\eta_{j}^{2}{C_{2}}^{2}} + {\zeta_{j}\eta_{j}2{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}} +} \end{matrix} \\ {{\xi_{j}\zeta_{j}2{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}} + {\gamma_{j}\xi_{j}\eta_{j}2{A_{2}}{C_{2}}\cos\;\varphi_{A_{2}C_{2}}}} \end{Bmatrix}}} & (5) \end{matrix}$ where coefficients A₁ and A₂ represent the amplitudes of the reference beams corresponding to the first and second frequency components of the input beam; coefficients B₁ and B₂ represent the amplitudes of background beams corresponding to reference beams A₁ and A_(2,) respectively; coefficients C₁ and C₂ represent the amplitudes of the return measurement beams corresponding to reference beams A₁ and A_(2,) respectively; P_(j) represents the integrated intensity of the first frequency component of the input beam in pulse j of the pulse sequence; and the values for ε_(j) and γ_(j) are listed in Table 1. The change in the values of ε_(j) and γ_(j) from 1 to −1 or from −1 to 1 correspond to changes in relative phases of respective reference and measurement beams. The coefficients ξ_(j), ζ_(j), and η_(j) represent effects of variations in properties of a conjugate set of four pinholes such as size and shape used in the generation of the spot on and/or in substrate 60 and the sensitivities of a conjugate set of four detector pixels corresponding to the spot on and/or in substrate 60 for the reference beam, the background beam, and the return measurement beam, respectively.

TABLE 1 j ε_(j) γ_(j) ε_(j)γ_(j) 1 1 1 1 2 −1 −1 1 3 −1 1 −1 4 1 −1 −1

It is assumed in Equation (5) that the ratio of |A₂|/|A₁| is not dependent on j or on the value of P_(j). In order to simplify the representation of S_(j) so as to project the important features without departing from either the scope or spirit of the present invention, it is also assumed in Equation (5) that the ratio of the amplitudes of the return measurement beams corresponding to A₂ and A₁ is not dependent on j or on the value of P_(j). However, the ratio |C₂|/|C₁| will be different from the ratio |A₂|/|A₁| when the ratio of the amplitudes of the measurement beam components corresponding to A₂ and A₁ are different from the ratio |A₂|/|A₁|.

Noting that cos Φ_(A) ₂ _(C) ₂ =±sin Φ_(A) ₁ _(C) ₁ by the control of the relative phase shifts between corresponding reference and return measurement beam components in beam 32, Equation (5) may be rewritten as

$\begin{matrix} {{S_{j} = {P_{j}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{j}^{2}\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right)} + {\zeta_{j}^{2}\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right)} + {\eta_{j}^{2}\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right)} +} \\ {{2\xi_{j}{\zeta_{j}\left( {{{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}} + {{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}}} \right)}} +} \end{matrix} \\ {{2\xi_{j}{\eta_{j}\begin{bmatrix} {{ɛ_{j}{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \\ {{\gamma_{j}\left( \frac{A_{2}}{A_{1}} \right)}\left( \frac{C_{2}}{C_{1}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} \end{bmatrix}}} +} \end{matrix} \\ {2\zeta_{j}{\eta_{j}\left\lbrack {{ɛ_{j}{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}} + {\gamma_{j}{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}}} \right\rbrack}} \end{Bmatrix}}},} & (6) \end{matrix}$ where the relationship cos Φ_(A) ₂ _(C) ₂ =sin Φ_(A) ₁ _(C) ₁ has been used without departing from either the scope or spirit of the present invention.

The change in phase Φ_(A) ₁ _(B) ₁ _(ε) _(j) Φ_(A) ₂ _(B) ₂ _(γ) _(j) for a change in ε_(j) and the change in phase Φ_(A) ₁ _(B) ₁ _(ε) _(j) Φ_(A) ₂ _(B) ₂ _(γ) _(j) for a change in γ_(j) may be different from π in embodiments depending on where and how the background beam is generated. It may be of value in evaluating the effects of the background beams to note that the factor cos Φ_(B) ₁ _(C) ₁ _(ε) _(j) may be written as cos [Φ_(A) ₁ _(C) ₁ +(Φ_(B) ₁ _(C) ₁ _(ε) _(j)−Φ_(A) ₁ _(C) ₁ )] where the phase difference (Φ_(B) ₁ _(C) ₁ _(ε) _(j) −Φ_(A) ₁ _(C) ₁ ) is the same as the phase φ_(A) ₁ _(B) ₁ _(ε) _(j) , i.e., cos Φ_(B) ₁ _(C) ₁ _(ε) _(j) =cos (Φ_(A) ₁ _(C) ₁ +Φ_(A) ₁ _(B) ₁ _(ε) _(j) ).

It is evident from inspection of Equation (6) that the term in Equation (6) corresponding to the component of conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ is a rectangular function that has a mean value of zero and is symmetric about j=2.5 since εj is symmetric about j=2.5. In addition the term in Equation (6) corresponding to the component of conjugated quadratures |C₁| sin Φ_(A) ₁ _(C) ₁ in Equation (6) is a rectangular function that has a mean value of zero and is antisymmetric about j=2.5 since γ_(j) is a antisymmetric function about j=2.5. Another important property by the design of the bi-homodyne detection method is that the conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ terms are orthogonal over the range of j=1,2,3,4 since ε_(j) and γ_(j) are orthogonal over the range of j=1,2,3,4,

${i.e.}\;,{{\sum\limits_{j = 1}^{4}{ɛ_{j}\gamma_{j}}} = 0.}$

Information about conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ is obtained using the symmetric and antisymmetric properties and orthogonality property of the conjugated quadratures terms in Equation (6) as represented by the following digital filters applied to the signals S_(j):

$\begin{matrix} {{F_{1}(S)} = {{\sum\limits_{j = 1}^{4}\;{ɛ_{j}\frac{S_{j}}{P_{j}^{\prime}\xi_{j}^{\prime 2}}}} = {{\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}\;{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}}} & (7) \\ {{{\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}\;{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{79mu}} & \; \\ {\mspace{85mu}{{\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}\;{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\eta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{79mu}} & \; \\ {\mspace{79mu}{{2{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{85mu}} & \; \\ {\mspace{76mu}{{2\left( \frac{A_{2}}{A_{1}} \right)\left( \frac{C_{2}}{C_{1}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}}} +}} & \; \\ {\mspace{70mu}{{2{A_{1}}{B_{1}}{\sum\limits_{j = 1}^{4}\;{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}}}} +}\mspace{40mu}} & \; \\ {\mspace{70mu}{{2{A_{2}}{B_{2}}{\sum\limits_{j = 1}^{4}\;{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}}}} +}\mspace{40mu}} & \; \\ {\mspace{70mu}{{2{B_{1}}{C_{1}}{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}}}} +}\mspace{65mu}} & \; \\ {\mspace{65mu}{{2{B_{2}}{C_{2}}{\sum\limits_{j = 1}^{4}\;{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}}}},}\mspace{34mu}} & \; \\ {{and}\mspace{644mu}} & \; \\ {{F_{2}(S)} = {{\sum\limits_{j = 1}^{4}\;{\gamma_{j}\frac{S_{j}}{P_{j}^{\prime}\xi_{j}^{\prime 2}}}} = {{\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}}} & (8) \\ {{{\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{79mu}} & \; \\ {\mspace{85mu}{{\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right){\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\eta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{79mu}} & \; \\ {\mspace{79mu}{{2{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}}} +}\mspace{85mu}} & \; \\ {\mspace{70mu}{{2\left( \frac{A_{2}}{A_{1}} \right)\left( \frac{C_{2}}{C_{1}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}}} +}} & \; \\ {\mspace{70mu}{{2{A_{1}}{B_{1}}{\sum\limits_{j = 1}^{4}\;{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}}}} +}\mspace{40mu}} & \; \\ {\mspace{70mu}{{2{A_{2}}{B_{2}}{\sum\limits_{j = 1}^{4}\;{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}}}} +}\mspace{40mu}} & \; \\ {\mspace{65mu}{{2{B_{1}}{C_{1}}{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}}}} +}\mspace{65mu}} & \; \\ {\mspace{65mu}{2{B_{2}}{C_{2}}{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}}}}\mspace{85mu}} & \; \end{matrix}$ wherein ξ′_(j) and P′_(j) are values used in the digital filters to represent ξ_(j) and P_(j).

The parameter

$\begin{matrix} \left\lbrack {\left( \frac{A_{2}}{A_{1}} \right)\left( \frac{C_{2}}{C_{1}} \right)} \right\rbrack & (9) \end{matrix}$ in Equations (7) and (8) needs to be determined in order complete the determination of a conjugated quadratures. The parameter given in Equation (9) can be measured for example by introducing π/2 phase shifts into the relative phase of the reference beam and the measurement beam and repeating the measurement for the conjugated quadratures. The ratio of the amplitudes of the conjugated quadratures corresponding to (sin Φ_(A) ₁ _(C) ₁ /cos Φ_(A) ₁ _(C) ₁ ) from the first measurement divided by the ratio of the amplitudes of the conjugated quadratures corresponding to (sin Φ_(A) ₁ _(C) ₁ /cos Φ_(A) ₁ _(C) ₁ ) from the second measurement is equal to

$\begin{matrix} {\left\lbrack {\left( \frac{A_{2}}{A_{1}} \right)\left( \frac{C_{2}}{C_{1}} \right)} \right\rbrack^{2}.} & (10) \end{matrix}$ Note that certain of the factors in Equations (7) and (8) have nominal values of 4 within scale factors, e.g.,

$\begin{matrix} {{{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}} \simeq 4},{{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}} \simeq 4.}} & (11) \end{matrix}$ The scale factors correspond to the average values for the ratios of ξ′_(j)/η_(j) and ξ′_(j)/ζ_(j), respectively, assuming that the average value of P_(j)/P′_(j)≅1. Certain other of the factors in Equations (7) and (8) have nominal values of zero, e.g.,

$\begin{matrix} \begin{matrix} {{{\sum\limits_{j = 1}^{4}{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} & {{{\sum\limits_{j = 1}^{4}{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} \\ {{{\sum\limits_{j = 1}^{4}{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\eta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} & \; \\ {{{\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} & {{{\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} \\ {{{\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\eta_{j}^{2}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0},} & \; \\ {{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\eta_{j}}{\xi_{j}^{\prime 2}} \right)}} \simeq 0.} & \; \end{matrix} & (12) \end{matrix}$ The remaining factors,

$\begin{matrix} \begin{matrix} {{\sum\limits_{j = 1}^{4}{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}}},{\sum\limits_{j = 1}^{4}{{ɛ_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}}},} \\ {{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}}},{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}}},} \\ {{\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{A_{1}B_{1}ɛ_{j}}}},{\sum\limits_{j = 1}^{4}{{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\xi_{j}\zeta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{A_{2}B_{2}\gamma_{j}}}},} \\ {{\sum\limits_{j = 1}^{4}{ɛ_{j}{\gamma_{j}\left( \frac{P_{j}}{P_{j}^{\prime}} \right)}\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{B_{1}C_{1}ɛ_{j}}}},{\sum\limits_{j = 1}^{4}{\left( \frac{P_{j}}{P_{j}^{\prime}} \right)\left( \frac{\zeta_{j}\eta_{j}}{\xi_{j}^{\prime\mspace{11mu} 2}} \right)\cos\;\varphi_{B_{2}C_{2}\gamma_{j}}}},} \end{matrix} & (13) \end{matrix}$ will have nominal magnitudes ranging from approximately zero to approximately 4 times a cosine factor and either the average value of factor (P_(j)/P′_(J))(ξ_(j)ζ_(j)/ξ′_(j) ²) or (P_(j)/P′_(J))(ζ_(j)η_(j)/ξ′_(j) ²) depending on the properties respective phases. For the portion of the back ground with phases that do not track to a first approximation the phases of the measurement beams, the magnitudes of all of the terms listed in the Equation (13) will be approximately zero. For the portion of the background with phases that do track to a first approximation the phases of the respective measurement beams, the magnitudes of the terms listed in Equation (13) will be approximately 4 times a cosine factor and either the average value of factor (P_(j)/P′_(J))(ξ_(j)ζ_(j)/ξ′_(j) ²) and or factor (P_(j)/P′_(J)) (ζ_(j)ηn_(j)/ξ′_(j) ²).

The two largest terms in Equations (7) and (8) are generally the terms that have the factors (|A₁|²+|A₂|²) and (|B₁|²+|B₂|²). However, the corresponding terms are substantially eliminated by selection of ξ′_(j) values for the terms that have (|A₁|²+|A₂|²) as a factor and by the design of ζ_(j) values for the terms that have (|B₁|²+|B₂|²) as a factor as shown in Equation (12).

The largest contribution from effects of background is represented by the contribution to the interference term between the reference beam and the portion of the background beam generated by the measurement beam component of beam 28. This portion of the effect of the background can be measured by measuring the corresponding conjugated quadratures of the portion of the background with the return measurement beam component of beam 32 set equal to zero, i.e., measuring the respective electrical interference signals S_(j) with substrate 60 removed and with either |A₂|=0 or |A₁|=0 and visa versa. The measured conjugated quadratures of the portion of the effect of the background can than used to compensate for the respective background effects beneficially in an end use application if required.

Information about the largest contribution from effects of background amplitude 2ξ_(j)ζ_(j)|A₁| |B₁| and phase Φ_(A) ₁ _(B) ₁ _(ε) _(j) , i.e., the interference term between the reference beam and the portion of background beam generated by the measurement beam component of beam 28, may be obtained by measuring S_(j) for j=1,2,3,4 as a function of relative phase shift between reference beam and the measurement beam component of beam 28 with substrate 60 removed and either |A₂|=0 or |A₁|=0 and visa versa and Fourier analyzing the measured values of S_(j). Such information can be used to help identify the origin of the respective background.

Other techniques may be incorporated into other embodiments to reduce and/or compensate for the effects of background beams without departing from either the scope or spirit of the present invention such as described in commonly owned U.S. Pat. No. 5,760,901 entitled “Method And Apparatus For Confocal Interference Microscopy With Background Amplitude Reduction and Compensation,” U.S. Pat. No. 5,915,048 entitled “Method and Apparatus for Discrimination In-Focus Images from Out-of-Focus Light Signals from Background and Foreground Light Sources,” and U.S. Pat. No. 6,480,285 B1 wherein each of the three patents are by Henry A. Hill. The contents of each of the three cited patents are herein incorporated in their entirety by reference.

The selection of values for ξ′_(j) based on information about coefficients ξ_(j) for j=1,2,3,4 that may be obtained by measuring the S_(j) for j=1,2,3,4 with only the reference beam present in the interferometer system. In certain embodiments, this may correspond simply blocking the measurement beam components of input beam 24 and in certain other embodiments, this may correspond to simply measuring the S_(j) for j=1,2,3,4 with substrate 60 removed. A test of the correctness of a set of values for ξ′_(j) is the degree to which the (|A₁|²+|A₂|²) terms in Equations (7) and (8) are zero.

Information about coefficients ξ_(j)η_(j) for j=1,2,3,4 may be obtained by scanning an artifact past the spots corresponding to the respective four conjugate detector pixels with either |A₂|=0 or |A₁|=0 and measuring the conjugated quadratures component 2|A₁| |C₁| cos Φ_(A) ₁ _(C) ₁ or 2|A₁| |C₁| sin Φ_(A) ₁ _(C) ₁ , respectively. A change in the amplitude of the 2|A₁| |C₁| cos Φ_(A) ₁ _(C) _(1 l) or 2|A₁| |C₁| sin Φ_(A) ₁ _(C) ₁ term corresponds to a variation in ξ_(j)η_(j) as a function of j. Information about the coefficients ξ_(j)η_(j) for j=1,2,3,4 may be used for example to monitor the stability of one or more elements of interferometer system 10.

The bi-homodyne detection method is a robust technique for the determination of conjugated quadratures of fields. First, the conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ are the primary terms in the digitally filtered values F₁ (S) and F₂ (S), respectively, since as noted in the discussion with respect to Equation (12), the terms with the factors (|A₁|²+|A₂|²) and (|B₁|²+|B₂|²) are substantially zero.

Secondly, the coefficients of |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₂| sin Φ_(A) ₁ _(C) ₁ terms in Equations (7) and (8) are identical. Thus highly accurate measurements of the interference terms between the return measurement beam and the reference beam with respect to amplitudes and phases, i.e., highly accurate measurements of conjugated quadratures of fields can be measured wherein first order variations in ξ_(j) and first order errors in normalizations such as (P_(j)/P′_(j)) and (ξ_(j) ²/ξ′_(j) ²) enter in only second or higher order. This property translates into a significant advantage. Also, the contributions to each component of the conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ from a respective set of four electrical interference signal values have the same window function and thus are obtained as jointly determined values.

Other distinguishing features of the bi-homodyne technique are evident in Equations (7) and (8): the coefficients of the conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ in Equations (7) and (8), respectively, corresponding to the first equation of Equations (11) are identical independent of errors in assumed values for ξ_(j) and η_(j); the coefficients of the conjugated quadratures |C₁| sin Φ_(A) ₁ _(C) ₁ and |C₁| cos Φ_(A) ₁ _(C) ₁ in Equations (7) and(8), respectively, corresponding to the fourth equation of Equations (12) are identical independent of errors in assumed values for ξ′_(j). Thus highly accurate values of the phases corresponding to conjugated quadratures can be measured with first order variations in ξ_(j) and first order errors in normalizations such as (P_(j)/P′_(j)) and (ξ_(j) ²/ξ′_(j) ²) enter in only through some high order effect.

It is also evident that since the conjugated quadratures of fields are obtained jointly when using the bi-homodyne detection method, there is a significant reduction in the potential for an error in tracking phase as a result of a phase redundancy unlike the situation possible in single-homodyne detection of conjugated quadratures of fields.

There are a number of advantages of the bi-homodyne detection as a consequence of the conjugated quadratures of fields being jointly acquired quantities. One advantage is a reduced sensitivity the effects of an overlay error of a spot in or on the substrate that is being imaged and a conjugate image of conjugate pixel of a multi-pixel detector during the acquisition of four electrical interference signal values of each spot in and/or on a substrate imaged using interferometric confocal microscopy. Overlay errors are errors in the set of four conjugate images of a respective set of conjugate detector pixels relative to the spot being imaged.

Another advantage is that when operating in the scanning mode there is a reduced sensitivity to effects of pinhole-to-pinhole variations in properties of a conjugate set of pinholes used in a confocal microscopy system that are conjugate to a spot in or on the substrate being imaged at different times during the scan.

Another advantage is that when operating in the scanning mode there is a reduced sensitivity to effects of pixel-to-pixel variation of properties within a set of conjugate pixels that are conjugate to a spot in or on the substrate being imaged at different times during the scan.

Another advantage is that when operating in the scanning mode there is reduced sensitivity to effects of pulse-to-pulse variations of a respective conjugate set of pulses of the input beam 24 to the interferometer system.

The pinholes and pixels of a multi-pixel detector of a set of conjugate pinholes and conjugate pixels of a multi-pixel detector may comprise contiguous pinholes of an array of pinholes and/or contiguous pixels of a multi-pixel detector or may comprise selected pinholes from an array of pinholes and/or pixels from an array of pixels wherein the separation between the selected pinholes is an integer number of pinhole separations and the separation between an array of respective pixels corresponds to an integer number of pixel separations without loss of lateral and/or longitudinal resolution and signal-to-noise ratios. The corresponding scan rate would be equal to the integer times the spacing of spots on the measurement object 60 conjugate to set of conjugate pinholes and/or set of conjugate pixels divided by the read out rate of the multi-pixel detector. This property permits a significant increase in through put for an interferometric confocal microscope with respect to the number of spots in and/or on a substrate imaged per unit time.

Referring to the quad-homodyne detection method, a set of four electrical interference signal values is obtained for each spot on and/or in substrate 60 being imaged with two pulse sequences from source 18 and beam-conditioner beam-conditioner 22. The set of four electrical interference signal values S_(j), j=1,2,3,4 used for obtaining conjugated quadratures of fields for a single a spot on and/or in a substrate being imaged is represented for the quad-homodyne detection within a scale factor by the formulae

$\begin{matrix} {{S_{1} = {P_{1}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{1}^{2}{A_{1}}^{2}} + {\zeta_{1}^{2}{B_{1}}^{2}} + {\eta_{1}^{2}{C_{1}}^{2}} + {\zeta_{1}\eta_{1}2{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{1}}} +} \\ {{\xi_{1}\zeta_{1}2{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{1}}} + {ɛ_{1}\xi_{1}\eta_{1}2{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \end{matrix} \\ {{\xi_{1}^{2}{A_{2}}^{2}} + {\zeta_{1}^{2}{B_{2}}^{2}} + {\eta_{1}^{2}{C_{2}}^{2}} + {\zeta_{1}\eta_{1}2{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{1}}} +} \end{matrix} \\ {{\xi_{1}\zeta_{1}2{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{1}}} + {\gamma_{1}\xi_{1}\eta_{1}2{A_{2}}{C_{2}}\cos\;\varphi_{A_{2}C_{2}}}} \end{Bmatrix}}},} & (14) \\ {{S_{2} = {P_{1}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{2}^{2}{A_{3}}^{2}} + {\zeta_{2}^{2}{B_{3}}^{2}} + {\eta_{2}^{2}{C_{3}}^{2}} + {\zeta_{2}\eta_{2}2{B_{3}}{C_{3}}\cos\;\varphi_{B_{3}C_{3}ɛ_{2}}} +} \\ {{\xi_{2}\zeta_{2}2{A_{3}}{B_{3}}\cos\;\varphi_{A_{3}B_{3}ɛ_{2}}} + {ɛ_{2}\xi_{2}\eta_{2}2{A_{3}}{C_{3}}\cos\;\varphi_{A_{3}C_{3}}} +} \end{matrix} \\ {{\xi_{2}^{2}{A_{4}}^{2}} + {\zeta_{2}^{2}{B_{4}}^{2}} + {\eta_{2}^{2}{C_{4}}^{2}} + {\zeta_{2}\eta_{2}2{B_{4}}{C_{4}}\cos\;\varphi_{B_{4}C_{4}\gamma_{2}}} +} \end{matrix} \\ {{\xi_{2}\zeta_{2}2{A_{4}}{B_{4}}\cos\;\varphi_{A_{4}B_{4}\gamma_{2}}} + {\gamma_{2}\xi_{2}\eta_{2}2{A_{4}}{C_{4}}\cos\;\varphi_{A_{4}C_{4}}}} \end{Bmatrix}}},} & (15) \\ {{S_{3} = {P_{2}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{1}^{2}{A_{1}}^{2}} + {\zeta_{1}^{2}{B_{1}}^{2}} + {\eta_{1}^{2}{C_{1}}^{2}} + {\zeta_{1}\eta_{1}2{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{3}}} +} \\ {{\xi_{1}\zeta_{1}2{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{3}}} + {ɛ_{3}\xi_{1}\eta_{1}2{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \end{matrix} \\ {{\xi_{1}^{2}{A_{2}}^{2}} + {\zeta_{1}^{2}{B_{2}}^{2}} + {\eta_{1}^{2}{C_{2}}^{2}} + {\zeta_{1}\eta_{1}2{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{3}}} +} \end{matrix} \\ {{\xi_{1}\zeta_{1}2{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{3}}} + {\gamma_{3}\xi_{1}\eta_{1}2{A_{2}}{C_{2}}\cos\;\varphi_{A_{2}C_{2}}}} \end{Bmatrix}}},} & (16) \\ {{S_{4} = {P_{2}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{2}^{2}{A_{3}}^{2}} + {\zeta_{2}^{2}{B_{3}}^{2}} + {\eta_{2}^{2}{C_{3}}^{2}} + {\zeta_{2}\eta_{2}2{B_{3}}{C_{3}}\cos\;\varphi_{B_{3}C_{3}ɛ_{4}}} +} \\ {{\xi_{2}\zeta_{2}2{A_{3}}{B_{3}}\cos\;\varphi_{A_{3}B_{3}ɛ_{4}}} + {ɛ_{4}\xi_{2}\eta_{2}2{A_{3}}{C_{3}}\cos\;\varphi_{A_{3}C_{3}}} +} \end{matrix} \\ {{\xi_{2}^{2}{A_{4}}^{2}} + {\zeta_{2}^{2}{B_{4}}^{2}} + {\eta_{2}^{2}{C_{4}}^{2}} + {\zeta_{2}\eta_{2}2{B_{4}}{C_{4}}\cos\;\varphi_{B_{4}C_{4}\gamma_{4}}} +} \end{matrix} \\ {{\xi_{2}\zeta_{2}2{A_{4}}{B_{4}}\cos\;\varphi_{A_{4}B_{4}\gamma_{4}}} + {\gamma_{4}\xi_{2}\eta_{2}2{A_{4}}{C_{4}}\cos\;\varphi_{A_{4}C_{4}}}} \end{Bmatrix}}},} & (17) \end{matrix}$ where coefficients A₁, A₂, A₃, and A₄ represent the amplitudes of the reference beams corresponding to the first, second, third, and fourth frequency components, respectively, of input beam 24; coefficients B₁, B₂, B₃, and B₄ represent the amplitudes of background beams corresponding to reference beams A₁, A₂, A₃, and A₄, respectively; coefficients C₁, C₂, C₃, and C₄ represent the amplitudes of the return measurement beams corresponding to reference beams A₁, A₂, A₃, and A₄, respectively; P₁ and P₂ represent the integrated intensities of the first frequency component in the first and second pulse sequences, respectively, of the input beam 24; and the values for ε_(j) and γ_(j) are listed in Table 1. The description of the coefficients ξ_(j), ζ_(j), and η_(j) for the quad-homodyne detection method is the same as the corresponding portion of the description given for ε_(j), ζ_(j), and η_(j) of the bi-homodyne detection method.

It is assumed in Equations (14), (15), (16), and (17) that the ratios of |A₂|/|A₁ and |A₄|/|A₃| are not dependent on j or the value of P_(j). In order to simplify the representation of S_(j) so as to project the important features without departing from either the scope or spirit of the present invention, it is also assumed in Equations (14), (15), (16), and (17) that the ratios of the amplitudes of the return measurement beams corresponding to |A₂|/|A₁| and |A₄|/|A₃| are not dependent on j or the value of P_(j). However, the ratios |C₂|/|C₁| and |C₄|/|C₃| will be different from the ratios |A₂|/|A₁| and |A₄|/|A₃|, respectively, when the ratio of the amplitudes of the measurement beam components corresponding to |A₂|/|A₁ and |A₄|/|A₃|, respectively, are different from the ratios |A₂ |/|A₁| and |A₄|/|A₃, respectively.

Noting that cos Φ_(A) ₂ _(C) ₂ =±sin Φ_(A) ₁ _(C) ₁ by the control of the relative phase shifts between corresponding reference and measurement beam components in beam 32, Equations (14), (15), (16), and (17) may be written, respectively, as

$\begin{matrix} {{S_{1} = {P_{1}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{1}^{2}\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right)} + {\zeta_{1}^{2}\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right)} + {\eta_{1}^{2}\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right)} +} \\ {{2\zeta_{1}{\eta_{1}\left\lbrack {{{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{1}}} + {{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{1}}}} \right\rbrack}} +} \end{matrix} \\ {{2\xi_{1}{\eta_{1}\begin{bmatrix} {{ɛ_{1}{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \\ {{\gamma_{1}\left( \frac{A_{2}}{A_{1}} \right)}\left( \frac{C_{2}}{C_{1}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} \end{bmatrix}}} +} \end{matrix} \\ {2\xi_{1}{\zeta_{1}\left\lbrack {{{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{1}}} + {{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{1}}}} \right\rbrack}} \end{Bmatrix}}},} & (18) \\ {{S_{2} = {P_{1}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{2}^{2}\left( {{A_{3}}^{2} + {A_{4}}^{2}} \right)} + {\zeta_{2}^{2}\left( {{B_{3}}^{2} + {B_{4}}^{2}} \right)} + {\eta_{2}^{2}\left( {{C_{3}}^{2} + {C_{4}}^{2}} \right)} +} \\ {{2\zeta_{2}{\eta_{2}\left\lbrack {{{B_{3}}{C_{3}}\cos\;\varphi_{B_{3}C_{3}ɛ_{2}}} + {{B_{4}}{C_{4}}\cos\;\varphi_{B_{4}C_{4}\gamma_{2}}}} \right\rbrack}} +} \end{matrix} \\ {{2\xi_{2}{\eta_{2}\left( \frac{A_{3}}{A_{1}} \right)}{\left( \frac{C_{3}}{C_{1}} \right)\begin{bmatrix} {{ɛ_{2}{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \\ {{\gamma_{2}\left( \frac{A_{4}}{A_{3}} \right)}\left( \frac{C_{4}}{C_{3}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} \end{bmatrix}}} +} \end{matrix} \\ {2\xi_{2}{\zeta_{2}\left\lbrack {{{A_{3}}{B_{3}}\cos\;\varphi_{A_{3}B_{3}ɛ_{2}}} + {{A_{4}}{B_{4}}\cos\;\varphi_{A_{4}B_{4}\gamma_{2}}}} \right\rbrack}} \end{Bmatrix}}},} & (19) \\ {{S_{3} = {P_{2}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{1}^{2}\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right)} + {\zeta_{1}^{2}\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right)} + {\eta_{1}^{2}\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right)} +} \\ {{2\zeta_{1}{\eta_{1}\left\lbrack {{{B_{1}}{C_{1}}\cos\;\varphi_{B_{1}C_{1}ɛ_{3}}} + {{B_{2}}{C_{2}}\cos\;\varphi_{B_{2}C_{2}\gamma_{3}}}} \right\rbrack}} +} \end{matrix} \\ {{2\xi_{1}{\eta_{1}\begin{bmatrix} {{ɛ_{3}{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \\ {{\gamma_{3}\left( \frac{A_{2}}{A_{1}} \right)}\left( \frac{C_{2}}{C_{1}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} \end{bmatrix}}} +} \end{matrix} \\ {2\xi_{1}{\zeta_{1}\left\lbrack {{{A_{1}}{B_{1}}\cos\;\varphi_{A_{1}B_{1}ɛ_{3}}} + {{A_{2}}{B_{2}}\cos\;\varphi_{A_{2}B_{2}\gamma_{3}}}} \right\rbrack}} \end{Bmatrix}}},} & (20) \\ {{S_{4} = {P_{2}\begin{Bmatrix} \begin{matrix} \begin{matrix} {{\xi_{2}^{2}\left( {{A_{3}}^{2} + {A_{4}}^{2}} \right)} + {\zeta_{2}^{2}\left( {{B_{3}}^{2} + {B_{4}}^{2}} \right)} + {\eta_{2}^{2}\left( {{C_{3}}^{2} + {C_{4}}^{2}} \right)} +} \\ {{2\zeta_{2}{\eta_{2}\left\lbrack {{{B_{3}}{C_{3}}\cos\;\varphi_{B_{3}C_{3}ɛ_{4}}} + {{B_{4}}{C_{4}}\cos\;\varphi_{B_{4}C_{4}\gamma_{4}}}} \right\rbrack}} +} \end{matrix} \\ {{2\xi_{2}{\eta_{2}\left( \frac{A_{3}}{A_{1}} \right)}{\left( \frac{C_{3}}{C_{1}} \right)\begin{bmatrix} {{ɛ_{4}{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +} \\ {{\gamma_{4}\left( \frac{A_{4}}{A_{3}} \right)}\left( \frac{C_{4}}{C_{3}} \right){A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} \end{bmatrix}}} +} \end{matrix} \\ {2\xi_{2}{\zeta_{2}\left\lbrack {{{A_{3}}{B_{3}}\cos\;\varphi_{A_{3}B_{3}ɛ_{4}}} + {{A_{4}}{B_{4}}\cos\;\varphi_{A_{4}B_{4}\gamma_{4}}}} \right\rbrack}} \end{Bmatrix}}},} & (21) \end{matrix}$ where the relationship cos Φ_(A) ₂ _(C) ₂ =sin Φ_(A) ₁ _(C) ₁ has been used without departing from either the scope or spirit of the present invention.

Information about the conjugated quadratures |C₁| cos Φ_(A) ₁ _(C) ₁ and |C₁| sin Φ_(A) ₁ _(C) ₁ are obtained using the symmetric and antisymmetric properties and orthogonality property of the conjugated quadratures as represented by the following digital filters applied to the signal values S_(j): j=1,2,3,4

$\begin{matrix} {{{F_{3}(S)} = {{\left( \frac{1}{P_{1}^{\prime}} \right)\left( {\frac{S_{1}}{\xi_{1}^{\prime 2}} - \frac{S_{2}}{\xi_{2}^{\prime 2}}} \right)} - {\left( \frac{1}{P_{2}^{\prime}} \right)\left( {\frac{S_{3}}{\xi_{1}^{\prime 2}} - \frac{S_{4}}{\xi_{2}^{\prime 2}}} \right)}}},} & (22) \\ {{F_{4}(S)} = {{\left( \frac{1}{P_{1}^{\prime}} \right)\left( {\frac{S_{1}}{\xi_{1}^{\prime 2}} - \frac{S_{2}}{\xi_{2}^{\prime 2}}} \right)} + {\left( \frac{1}{P_{2}^{\prime}} \right){\left( {\frac{S_{3}}{\xi_{1}^{\prime 2}} - \frac{S_{4}}{\xi_{2}^{\prime 2}}} \right).}}}} & (23) \end{matrix}$ The description of ξ′_(j) and P′_(j) for the quad-homodyne detection method is the same as the corresponding description given for ξ′_(j) and P′_(j) in the bi-homodyne detection method. Using Equations (18), (19), (20), (21), (22), and (23), the following expressions are obtained for the filtered quantities containing components of the conjugated quadratures

$\begin{matrix} {{C_{1}}\cos\;\varphi_{A_{1}C_{1}}\mspace{14mu}{and}} & \; \\ {{C_{1}}\;\sin\;\varphi_{A_{1}C_{1}}\text{:}} & \; \\ {{F_{3}(S)} = {{\left( {\frac{P_{1}}{P_{1}^{\prime}} - \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right)\left( \frac{\xi_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{A_{3}}^{2} + {A_{4}}^{2}} \right)\left( \frac{\xi_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & (24) \\ {{~~~~~~~~~~~~~~~~~}{{\left( {\frac{P_{1}}{P_{1}^{\prime}} - \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right)\left( \frac{\zeta_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{B_{3}}^{2} + {B_{4}}^{2}} \right)\left( \frac{\zeta_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & \; \\ {{~~~~~~~~~~~~~~~~~}{{\left( {\frac{P_{1}}{P_{1}^{\prime}} - \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right)\left( \frac{\eta_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{C_{3}}^{2} + {C_{4}}^{2}} \right)\left( \frac{\eta_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & \; \\ {{~~~~~~~~~~~~~~~~~}{2{\left( {\frac{P_{1}}{P_{1}^{\prime}} + \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {\left( \frac{\xi_{1}\eta_{1}}{\xi_{1}^{\prime 2}} \right) + {\left( \frac{\xi_{2}\eta_{2}}{\xi_{2}^{\prime 2}} \right)\left( \frac{A_{3}}{A_{1}} \right)\left( \frac{C_{3}}{C_{1}} \right)}} \right\rbrack}}} & \; \\ {{~~~~~~~~~~~~~~~~~}{{{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~}{2\left( {\frac{P_{1}}{P_{1}^{\prime}} - \frac{P_{2}}{P_{2}^{\prime}}} \right)\left( \frac{A_{2}}{A_{1}} \right){\left( \frac{C_{2}}{C_{1}} \right)\begin{bmatrix} {\left( \frac{\xi_{1}\eta_{1}}{\xi_{1}^{\prime 2}} \right) +} \\ {\left( \frac{\xi_{2}\eta_{2}}{\xi_{2}^{\prime 2}} \right)\left( \frac{A_{4}}{A_{2}} \right)\left( \frac{C_{4}}{C_{2}} \right)} \end{bmatrix}}}} & \; \\ {\mspace{85mu}{{{A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{1}B_{1}ɛ_{1}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{1}B_{1}ɛ_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{A_{1}}{B_{1}}} -}} & \; \\ {{~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{3}B_{3}ɛ_{2}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{3}B_{3}ɛ_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{A_{3}}{B_{3}}} +}} & \; \\ {\mspace{85mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{2}B_{2}\gamma_{1}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{2}B_{2}\gamma_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{A_{2}}{B_{2}}} -}} & \; \\ {\mspace{85mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{4}B_{4}\gamma_{2}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{4}B_{4}\gamma_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{A_{4}}{B_{4}}} +}} & \; \\ {\mspace{85mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{1}C_{1}ɛ_{1}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{1}C_{1}ɛ_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{B_{1}}{C_{1}}} -}} & \; \\ {\mspace{85mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{3}C_{3}ɛ_{2}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{3}C_{3}ɛ_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{B_{3}}{C_{3}}} +}} & \; \\ {\mspace{85mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{2}C_{2}\gamma_{1}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{2}C_{2}\gamma_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{B_{2}}{C_{2}}} -}} & \; \\ {{~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{4}C_{4}\gamma_{2}}} - {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{4}C_{4}\gamma_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{B_{4}}{C_{4}}},}} & \; \\ {{F_{4}(S)} = {{\left( {\frac{P_{1}}{P_{1}^{\prime}} + \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{A_{1}}^{2} + {A_{2}}^{2}} \right)\left( \frac{\xi_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{A_{3}}^{2} + {A_{4}}^{2}} \right)\left( \frac{\xi_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & (25) \\ {{~~~~~~~~~~~~~~~~~~~}{{\left( {\frac{P_{1}}{P_{1}^{\prime}} + \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{B_{1}}^{2} + {B_{2}}^{2}} \right)\left( \frac{\zeta_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{B_{3}}^{2} + {B_{4}}^{2}} \right)\left( \frac{\zeta_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{\left( {\frac{P_{1}}{P_{1}^{\prime}} + \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {{\left( {{C_{1}}^{2} + {C_{2}}^{2}} \right)\left( \frac{\eta_{1}^{2}}{\xi_{1}^{\prime 2}} \right)} - {\left( {{C_{3}}^{2} + {C_{4}}^{2}} \right)\left( \frac{\eta_{2}^{2}}{\xi_{2}^{\prime 2}} \right)}} \right\rbrack} +}} & \; \\ {\mspace{95mu}{2{\left( {\frac{P_{1}}{P_{1}^{\prime}} - \frac{P_{2}}{P_{2}^{\prime}}} \right)\left\lbrack {\left( \frac{\xi_{1}\eta_{1}}{\xi_{1}^{\prime 2}} \right) + {\left( \frac{\xi_{2}\eta_{2}}{\xi_{2}^{\prime 2}} \right)\left( \frac{A_{3}}{A_{1}} \right)\left( \frac{C_{3}}{C_{1}} \right)}} \right\rbrack}}} & \; \\ {{~~~~~~~~~~~~~~~~~~~~}{{{A_{1}}{C_{1}}\cos\;\varphi_{A_{1}C_{1}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{2\left( {\frac{P_{1}}{P_{1}^{\prime}} + \frac{P_{2}}{P_{2}^{\prime}}} \right)\left( \frac{A_{2}}{A_{1}} \right){\left( \frac{C_{2}}{C_{1}} \right)\begin{bmatrix} {\left( \frac{\xi_{1}\eta_{1}}{\xi_{1}^{\prime 2}} \right) +} \\ {\left( \frac{\xi_{2}\eta_{2}}{\xi_{2}^{\prime 2}} \right)\left( \frac{A_{4}}{A_{2}} \right)\left( \frac{C_{4}}{C_{2}} \right)} \end{bmatrix}}}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{{A_{1}}{C_{1}}\sin\;\varphi_{A_{1}C_{1}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{1}B_{1}ɛ_{1}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{1}B_{1}ɛ_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{A_{1}}{B_{1}}} -}} & \; \\ {\mspace{95mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{3}B_{3}ɛ_{2}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{3}B_{3}ɛ_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{A_{3}}{B_{3}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{2}B_{2}\gamma_{1}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{2}B_{2}\gamma_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{A_{2}}{B_{2}}} -}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{A_{4}B_{4}\gamma_{2}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{A_{4}B_{4}\gamma_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{A_{4}}{B_{4}}} +}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{1}C_{1}ɛ_{1}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{1}C_{1}ɛ_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{B_{1}}{C_{1}}} -}} & \; \\ {\mspace{95mu}{{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{3}C_{3}ɛ_{2}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{3}C_{3}ɛ_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{B_{3}}{C_{3}}} +}} & \; \\ {\mspace{95mu}{{2\;\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{2}C_{2}\gamma_{1}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{2}C_{2}\gamma_{3}}}} \right)\frac{\xi_{1}\zeta_{1}}{\xi_{1}^{\prime 2}}{B_{2}}{C_{2}}} -}} & \; \\ {{~~~~~~~~~~~~~~~~~~~}{2\left( {{\frac{P_{1}}{P_{1}^{\prime}}\cos\;\varphi_{B_{4}C_{4}\gamma_{2}}} + {\frac{P_{2}}{P_{2}^{\prime}}\cos\;\varphi_{B_{4}C_{4}\gamma_{4}}}} \right)\frac{\xi_{2}\zeta_{2}}{\xi_{2}^{\prime 2}}{B_{4}}{{C_{4}}.}}} & \; \end{matrix}$

The parameters

$\begin{matrix} {\left\lbrack {\left( \frac{A_{2}}{A_{1}} \right)\left( \frac{C_{2}}{C_{1}} \right)} \right\rbrack,} & (26) \\ {{\left( \frac{A_{4}}{A_{2}} \right)\left( \frac{C_{4}}{C_{2}} \right)},} & (27) \\ \left\lbrack {\left( \frac{A_{3}}{A_{1}} \right)\left( \frac{C_{3}}{C_{1}} \right)} \right\rbrack & (28) \end{matrix}$ need to be determined in order to complete the determination of a conjugated quadratures for certain end use applications. The parameters given by Equations (26), (27), and (28) can for example be measured by procedures analogous to the procedure described for the bi-homodyne detection method with respect to measuring the quantity specified by Equation (9).

The remaining description of the quad-homodyne detection method is the same as corresponding portion of the description given for the bi-homodyne detection method.

It is also evident that since the conjugated quadratures of fields are obtained jointly when using the quad-homodyne detection, there is a significant reduction in the potential for an error in tracking phase as a result of a phase redundancy unlike the situation possible in single-homodyne detection of conjugated quadratures of fields.

There are a number of advantages of the quad-homodyne detection as a consequence of the conjugated quadratures of fields being jointly acquired quantities.

One advantage is a reduced sensitivity the effects of an overlay error of a spot in or on the substrate that is being imaged and a conjugate image of a pixel of a conjugate set of pixels of a multi-pixel detector during the acquisition of the four electrical interference signal values of each spot in and/or on a substrate imaged using interferometric confocal microscopy. Overlay errors are errors in the set of four conjugate images of a respective set of conjugate detector pixels relative to the spot being imaged.

Another advantage is that when operating in the scanning mode there is reduced sensitivity to effects of pulse-to-pulse variations of a respective conjugate set of pulses of the input beam 24 to the interferometer system.

Another advantage is that when operating in the scanning mode there is an increase in throughput since only one pulse of the source is required to generate the at least four electrical interference values.

A first embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 of the first embodiment shown schematically in FIG. 2 a. Interferometer 10 comprises an interferometer such as described in commonly owned U.S. Provisional Patent Application No. 60/447,254 (ZI-40) entitled “Transverse Differential Interferometric Confocal Microscopy” and U.S patent application Ser. No. 10/778,371, filed Feb. 13, 2004 (ZI-40) also entitled “Transverse Differential Interferometric Confocal Microscopy” both of which are by Henry A. Hill. The contents of the U.S. Provisional Patent Application and the U.S. Patent Application are herein incorporated in their entirety by reference.

Interferometer 10 of the first embodiment comprises a first imaging system generally indicated as numeral 110, pinhole array beam-splitter 112, detector 70, and a second imaging system generally indicated as numeral 210. The second imaging system 210 is low power microscope having a large working distance, e.g. Nikon ELWD and SLWD objectives and Olympus LWD, ULWD, and ELWD objectives. The first imaging system 110 comprises an interferometric confocal microscopy system described in part in commonly owned U.S. Provisional Application No. 60/442,982 [ZI-45] entitled “Interferometric Confocal Microscopy Incorporating Pinhole Array Beam-Splitter” and U.S. patent application Ser. No. 10/765,229, filed Jan. 27, 2004 (ZI-45) and also entitled “Interferometric Confocal Microscopy Incorporating Pinhole Array Beam-Splitter” both of which are by Henry A. Hill. The contents of both of the U.S. Provisional Patent Application and the U.S. Patent Application are herein incorporated in their entirety by reference.

First imaging system 110 is shown schematically in FIG. 2 b. The imaging system 110 is a catadioptric system such as described in commonly owned U.S. Pat. No. 6,552,852 B2 (ZI-38) entitled “Catoptric and Catadioptric Imaging System” and commonly owned U.S. patent application Ser. No. 10/366,651 (ZI-43) entitled “Catoptric And Catadioptric Imaging Systems” wherein both of the patent applications are by Henry A. Hill, the contents of the two cited patent applications incorporated herein in their entirety by reference.

Catadioptric imaging system 110 comprises catadioptric elements 140 and 144, beam-splitter 148, and convex lens 150. Surfaces 142B and 146B are concave spherical surfaces with nominally the same radii of curvature and the conjugate points with respect to beam-splitter 148 of respective centers of curvature of surfaces 142B and 146B are the same. Surfaces 142A and 146A are convex spherical surfaces with nominally the same radii of curvature. The centers of curvature of surfaces 142A and 146A are transversely shifted in a plane parallel to the plane of beam-splitter 148 by a small radial displacements Δr₁ and Δr₂, respectively, and longitudinal displacements Δz₁ and Δz₂, respectively with respect to the centers of curvature of surfaces 146B and 142B, respectively. The relative orientation of the displacement vectors corresponding to Δr₁ and Δr₂ is chosen to minimize the affects of spurious beams as subsequently described. The relative displacements Δz₁ and Δz₂ are selected to optimize the performance of interferometer 10 with respect to acquisition of information about the surface of substrate 60.

As a result of the small displacements just mentioned, the conjugate of the center of curvature of surface 142A, as seen through beam splitter 148, does not coincide with the center of curvature of surface 146A. (Or, equivalently, the conjugate of the center of curvature of surface 146A, as seen through beam splitter 148, does not coincide with the center of curvature of surface 142A.) Rather, those two points are displaced by an amount determined by the small displacement of the two surfaces 142A and 146A relative to each other. And the vector of that displacement has a component that is normal to the plane of beam splitter 148 and a component that is parallel to the plane of beam splitter 148. Stated a different way, the line connecting the centers of curvature of the two surfaces 142A and 146A is not normal to the plane of beam splitter 148 but rather it diverges from the normal by a small angle that is determined by the relative displacement of the two surfaces.

The center of curvature of convex lens 150 is the same as the center of curvature of surfaces 142B. The radius of curvature of surface 146B is selected so as to minimize the loss in usable solid angle of the imaging system 110 and to produce a working distance for imaging system 110 acceptable for an end use application, e.g., of the order of a mm. The radius of curvature of convex lens 150 is selected so that off-axis aberrations of the catadioptric imaging system 110 are compensated. The medium of elements 140 and 144 may be for example CaF₂, fused silica or commercially available glass such as SF11. The medium of convex lens 150 may be for example CaF₂, fused silica, YAG, or commercially available glass such as SF11. An important consideration in the selection of the media of elements 140 and 144 and convex lens 150 will the transmission properties for the frequencies of beam 24.

Convex lens 152 has a center of curvature the same as the center of curvature of convex lens 150. Convex lenses 150 and 152 are bonded together with pinhole beam-splitter 112 in between. Pinhole array beam-splitter 112 is shown in FIG. 2 c. The pattern of pinholes in pinhole array beam-splitter is chosen to match the requirements of an end use application. An example of a pattern is a two dimensional array of equally spaced pinholes in two orthogonal directions. The pinholes may comprise circular apertures, rectangular apertures, or combinations thereof such as described in commonly owned U.S. patent application Ser. No. 09/917,402 (ZI-15) entitled “Multiple-Source Arrays for Confocal and Near-field Microscopy” by Henry A. Hill and Kyle Ferrio of which the contents thereof are incorporated herein in their entirety by reference. The spacing between pinholes of pinhole array beam-splitter 112 is shown in FIG. 2 c as b with aperture size a.

Input beam 24 is reflected by mirror 54 to pinhole beam-splitter 112 where a first portion thereof is transmitted as reference beam components of output beam components 130A and 130B (see FIG. 2 a) and a second portion thereof scattered as measurement beam components of beam components 126A and 126B. The measurement beam components of beam components 126A and 126B are imaged as measurement beam components of beam components 128A and 128B to an array of image spots in image planes displaced from the surface of substrate 60.

The arrays of image spots in the image planes displaced from the surface of substrate 60 comprises a first and second array of image spots with the second array of image spots transversely and longitudinally displaced with respect to the first array of image spots. A corresponding pair of spots 164 and 166 of the first and second arrays of image spots are shown diagrammatically in FIGS. 2 d and 2 e for the case where the displacements of convex surfaces 142A and 146A lie in the x-z plane of FIGS. 2 d and 2 e and (−2x₁,0,2z₁) and (−2x₂,0,−2z₂) are the respective locations of the pair of image spots. The displacements of surfaces 142A and 146A are restricted to the x-z plane in order to simplify the description of the important aspects without limiting the scope or the spirit of the invention. The displacement of center of curvature of surface 142A is in the negative x direction and positive z direction. The displacement of center of curvature of surface 146A is in the positive x direction and negative z direction. The displacements of the pair of image spots 164 and 166 are a consequence of the displacements of centers of curvature of surfaces 142A and 146A. An example of a path of a beam contributing to image spot 164 is beam 126E and an example of a path of a beam contributing to image spot 166 is beam 126F (see FIG. 2 d). Note that in FIG. 2 d, the left side of the catadioptric element illustrates how the reflected measurement beam (i.e., reflected off of beam splitter 148) is focused onto spot 164 and the right side of the catadioptric element illustrates how the transmitted measurement beam is focused onto spot 166. There is also a reflected measurement beam on the right side which, like the reflected measurement beam on the left side, is focused onto spot 164. Similarly, there is a transmitted measurement beam on the left side, which, like the transmitted measurement beam on the right side, is focused onto spot 166.

A portion of beam 126E is also reflected twice by beam-splitter 148 and once by convex surface 142A to form an image spot 184 (see FIG. 2 d) at location (−2x₁/n,0,−2z₁/n) where n is the index of refraction of convex lenses 150 and 152. In addition, a portion of beam 126F is transmitted twice by beam-splitter 148 and reflected once by convex surface 146A to form an image spot 186 (see FIG. 2 d) at location (2x₂/n,0,2z₂/n).

Next consider the affects of catadioptric imaging system 110 on the portions of the beams comprising image spots 164 and 166 that are reflected by the surface of substrate 60. The reflected image spots are represented by image sources 1641 and 166I as shown in the figure. The reflected portions are part of the return measurement beam components of beams 128A and 128B and imaged by catadioptric imaging system 110 to four spots 190, 192, 194, and 196 in the space of pinhole array beam-splitter 112 (see FIG. 2 e). The locations of the respective spots with respect to the pinhole source of beams 126E and 126F are (0,0, [2h₁−4z₁]/2), (0,0, [2h₂+4z₁]/2), (−2[x₁+x₂]/n,0,2h₂/n), and (2[x₁+x₂]/n,0,2h₁/n) where h₁ and h₂ are the heights of image spots 164 and 166, respectively, below the surface of substrate 60. The surface of substrate 60 is not indicated in FIG. 2 e. The intersection of the vertical and horizontal dot-dash-dot lines located between image spots 164 and 166 in FIG. 2 e corresponds to center of curvature of concave surface 146B (see FIG. 2 b).

Portions of the beams forming image spot 190 and 192 are transmitted by the pinhole corresponding to pinhole source of beams 126E and 126F as a component of beam components 130A and 130B. The beams forming image spots 184, 186, 194, and 196 are not transmitted by the pinhole corresponding to the pinhole source of beams 126E and 126F as a property of confocal imaging system 110. To further reduce the effects of spurious beams generated by reflection of portions of the beams forming image spots 184, 186, 194, and 196 by pinhole array 112, the respective surface of pinhole array 112 is coated with an anti-reflective coating.

The effects of spurious beams generated by reflection of portions of the beams forming image spots 184, 186, 194, and 196 by pinhole array 112 are reduced when the displacements in the x-y plane of the centers of curvature of convex surfaces 142A and 146A are selected to be orthogonal, e.g., the displacements of the centers of convex surfaces 142A and 146A are in the x-z and y-z planes, respectively. The corresponding displacements are Δx₁ and Δy₂. The resulting locations of image spots 184, 186, 194, and 196 in x-y plane are shown schematically in FIG. 2 f. With the image spot pattern shown in FIG. 2 f, only the reflection by the low reflectivity surface of pinhole array 112 of beams forming image spots 194 and 196 will generate secondary spurious beams wherein portions thereof can be subsequently transmitted by the pinhole corresponding to the pinhole source of beams 126E and 126F. The amplitudes of the effects of the secondary spurious beams will be less than or of the order of 2% of the amplitudes of the portions of the beams forming spots 190 and 192 transmitted by the corresponding pinhole in pinhole array 112. For the resulting locations of image spots 184, 186, 192, 194, and 196 shown schematically in FIG. 2 c, the amplitudes of the effects of the corresponding spurious beams will be less than or of the order of 40% of the amplitude of the portions of the beams forming spots 190 and 192 transmitted by the corresponding pinhole in pinhole array 112.

The description of the affects of catadioptric imaging system 110 on portions of the beams comprising image spots 164 and 166 that are scattered by sub-wavelength artifacts and/or defects on the surface of substrate 60 is based on an analysis that is a variant of the analysis forming the basis of the description given of the affects of catadioptric imaging system 110 on the portions of the beams comprising image spots 164 and 166 that are reflected by the surface of substrate 60.

The next step is the imaging of output beam components 130A and 130B by imaging system 210 to an array of spots that coincide with the pixels of a multi-pixel detector 70 such as a CCD to generate an array of electrical interference signals 72. The array of electrical interference signals is transmitted to signal processor and controller 80 for subsequent processing for an array of conjugated quadratures.

The description of input beam 24 is the same as corresponding portions of the description given for input beam 24 of FIG. 1 a with beam-conditioner 22 configured as a two-frequency generator and frequency-shifter shown in FIG. 1 c. Input beam 24 comprises two components that have different frequencies and have the same state of plane polarization. The frequency of each component of input beam 24 is shifted between different frequency values by beam-conditioner 22 according to control signal 74 generated by electronic processor and controller 80. Beam 20 comprises a single frequency component.

The conjugated quadratures of fields of the return measurement beams are obtained using the bi-homodyne detection method wherein sets of four measurements of the electrical interference signals 72 are made. An array of conjugated quadratures of fields is measured interferometrically by interferometer confocal imaging system 10 wherein each conjugated quadratures comprises a difference of conjugated quadratures of fields of beams scattered/reflected from a pair of spots in or on a substrate. The array of conjugated quadratures is measured jointly, i.e., simultaneously, and the components of each conjugated quadratures are measured jointly.

The relative phases of the beams subsequently scattered/reflected by the pair of spots in or on a substrate are adjusted by the control of a single interferometer system parameter, the relative optical path lengths of measurement and return measurement beam components respective spots of the pairs of spots in or on substrate 60. The relative phases are adjusted by making a longitudinal displacement of convex surface 142A with respect to convex surface 146A. This may be achieved by changing the thickness of element 140 relative to the thickness of element 144. Alternatively, a thin layer may be added to surface 142A and/or surface 146A. Another example of a technique to introduce a change in relative phase is to add a concave reflecting surface next to convex surface 142A with an air gap such that the air gap thickness may be adjusted. In the latter example, the convex surface 142A would be anti-reflective coated.

The measured conjugate quadratures in the first embodiment are proportional to components of complex amplitude V₂ (h₁,z₁,h₂, z₂,χ) can to a good approximation be written as

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},\chi} \right)} = {{R_{1}^{1/2}{j_{0}\left( {\alpha\left( {h_{1} - {2z_{1}}} \right)} \right\rbrack}{\mathbb{e}}^{{- {\mathbb{i}}}\;{\beta{({h_{1} - {2z_{1}}})}}}} +}} \\ {{\mathbb{e}}^{{\mathbb{i}}\;\chi}R_{2}^{1/2}{j_{0}\left\lbrack {\alpha\left( {h_{2} + {2z_{2}}} \right)} \right\rbrack}{\mathbb{e}}^{{- {\mathbb{i}}}\;{\beta{({h_{2} + {2z_{2}}})}}}} \end{matrix} & (29) \end{matrix}$ where R₁ ^(1/2) and R₂ ^(1/2) are the complex reflectivity coefficients of the surface of substrate 60 for the fields of the beams forming spots 190 and 192, respectively, α=2k(1−cos θ₀),  (30) β=2k(1+cos θ₀),  (31) j_(p) (x) is the spherical Bessel function of order p=0,1, 2, . . . , sin θ₀ is the numerical aperture of the imaging system in the image space of substrate 60, and χ is the relative phase of the components of beams forming spots 190 and 192 at a corresponding pinhole of pinhole array beam-splitter 112. Descriptions of derivations that form the basis for the two terms in Equation (29) may be found in references such as the book edited by T. Wilson, Confocal Microscopy, Academic Press (1990), the contents of which are herein incorporated in their entirety by reference.

The coefficient exp(iχ) in Equation (29) is written as a coefficient with a magnitude of 1 to a good approximation: the pinhole sources for each of the beams forming spots 164 and 166 (see FIG. 2 d) are the same pinhole of pinhole array beam-splitter 112, the pinholes performing the confocal spatial filtering of the beams forming spots 190 and 192 are the same pinhole of pinhole array beam-splitter 112 (see FIGS. 2 c and 2 e), and the portions of the beams forming spots 190 and 192 subsequently spatially filtered, i.e., transmitted, by the same pinhole are detected by the same pixel of detector 70. This feature represents one of the distinguishing advantages of the imaging system 110.

The second embodiment comprises the first embodiment configured for operation in a dark field mode. In the dark field configuration of the second embodiment, χ=π.  (32) and Equation (29) assumes the form

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},{\chi = \pi}} \right)} = {{R_{1}^{1/2}{j_{0}\left( {\alpha\left( {h_{1} - {2z_{1}}} \right)} \right\rbrack}{\mathbb{e}}^{{- {\mathbb{i}}}\;{\beta{({h_{1} - {2z_{1}}})}}}} -}} \\ {R_{2}^{1/2}{j_{0}\left\lbrack {\alpha\left( {h_{2} + {2z_{2}}} \right)} \right\rbrack}{{\mathbb{e}}^{{- {\mathbb{i}}}\;{\beta{({h_{2} + {2z_{2}}})}}}.}} \end{matrix} & (33) \end{matrix}$ The remainder of the description of the properties of the measured conjugated quadratures will be with respect to Equation (33) and the second embodiment. The extension of the description to a description of the measured conjugated quadratures of the first embodiment and Equation (29) will be apparent to one skilled in the art.

Properties of Equation (33) are easily recognized by examination of low order terms of h₁ and h₂ in a power series representation, i.e.,

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},{\chi = \pi}} \right)} = \left\lbrack {{R_{1}^{1/2}{j_{0}\left( {{- 2}\alpha\; z_{1}} \right)}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} -} \right.} \\ {\left. {R_{2}^{1/2}{j_{0}\left( {2\alpha\; z_{2}} \right)}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}} \right\rbrack -} \\ {\alpha\left\lbrack {{R_{1}^{1/2}{j_{1}\left( {{- 2}\alpha\; z_{1}} \right)}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}h_{1}} -} \right.} \\ {\left. {R_{2}^{1/2}{j_{1}\left( {2\alpha\; z_{2}} \right)}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}h_{2}} \right\rbrack -} \\ {{\mathbb{i}}\;{\beta\left\lbrack {{R_{1}^{1/2}{j_{0}\left( {{- 2}\alpha\; z_{1}} \right)}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}h_{1}} -} \right.}} \\ {\left. {R_{2}^{1/2}{j_{0}\left( {2\alpha\; z_{2}} \right)}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}h_{2}} \right\rbrack.} \end{matrix} & (34) \end{matrix}$ For the non-limiting assumption z₁=z₂, Equation (34) simplifies to the expression

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},{\chi = \pi}} \right)} = {{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left\lbrack {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} -} \right.}} \\ {\left. {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}} \right\rbrack +} \\ {\alpha\;{{j_{1}\left( {2\;\alpha\; z_{1}} \right)}\left\lbrack {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}h_{1}} +} \right.}} \\ {\left. {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}h_{2}} \right\rbrack -} \\ {{\mathbb{i}}\;\beta\;{{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left\lbrack {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}h_{1}} -} \right.}} \\ {\left. {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}h_{2}} \right\rbrack.} \end{matrix} & (35) \end{matrix}$ By combining terms of h₁ and h₂, Equation (35) reduces to

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},{\chi = \pi}} \right)} = {{{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left\lbrack {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} - {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{2}}}} \right\rbrack} +}} \\ {{R_{1}^{1/2}{{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}\left\lbrack {{\alpha\;{j_{1}\left( {2\;\alpha\; z_{1}} \right)}} - {{\mathbb{i}}\;\beta\;{j_{0}\left( {2\;\alpha\; z_{1}} \right)}}} \right\rbrack}h_{1}} +} \\ {R_{2}^{1/2}{{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}\left\lbrack {{\alpha\;{j_{1}\left( {2\;\alpha\; z_{1}} \right)}} + {{\mathbb{i}}\;\beta\;{j_{0}\left( {2\;\alpha\; z_{1}} \right)}}} \right\rbrack}{h_{2}.}} \end{matrix} & (36) \end{matrix}$

Another useful form of Equation (35) is obtained by writing it in terms of (h₁+h₂)/2 and (h₁−h₂)/2.

$\begin{matrix} \begin{matrix} {{V_{2}\left( {h_{1},z_{1},h_{2},z_{2},{\chi = \pi}} \right)} = {{{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left\lbrack {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} - {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}}} \right\rbrack} +}} \\ {{\begin{Bmatrix} {{\alpha\;{j_{1}\left( {2\;\alpha\; z_{1}} \right)}\left( {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} + {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}}} \right)} -} \\ {{\mathbb{i}}\;\beta\;{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left( {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} - {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}}} \right)} \end{Bmatrix}\left( \frac{h_{1} + h_{2}}{2} \right)} +} \\ {\frac{1}{2}\begin{Bmatrix} {{\alpha\;{j_{1}\left( {2\;\alpha\; z_{1}} \right)}\left( {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} - {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}}} \right)} -} \\ {{\mathbb{i}}\;\beta\;{j_{0}\left( {2\;\alpha\; z_{1}} \right)}\left( {{R_{1}^{1/2}{\mathbb{e}}^{{\mathbb{i}}\; 2\;\beta\; z_{1}}} + {R_{2}^{1/2}{\mathbb{e}}^{{- {\mathbb{i}}}\; 2\;\beta\; z_{1}}}} \right)} \end{Bmatrix}{\left( {h_{1} - h_{2}} \right).}} \end{matrix} & (37) \end{matrix}$

One important application of the second embodiment is the determination of differences in height of a first region comprising a feature or an artifact relative to a neighboring second region used as a reference region. The first and second regions would have heights corresponding to h₁ and h₂, respectively. The factors (h₁+h₂)/2 and (h₁−h₂) represent the average height and the difference in heights of the two regions on the surface of substrate 60 corresponding to beams forming spots 190 and 192. The coefficients of h₁ and h₂ in Equation (36) and coefficients of (h₁+h₂)/2 and (h₁−h₂) in Equation (37) can be independently measured by measuring V₂ (h₁,z₁, h₂, z₂,χ=π) with substrate 60 stationary in the x-y plane and with transducers 86A and 86B introducing a scan in z to obtain measurements at at least two different positions (thereby changing h₁ and h₂ and causing the average (h₁+h₂)/2 to change but keeping the difference (h₁−h₂) the same) and scans in orientation of substrate 60 for known values of R₁ ^(1/2) and R₂ ^(1/2) (wherein the scans are achieved by rotating the substrate to obtain at least two measurements in which the average (h₁+h₂)/2 does not change but the slope (h₁−h₂) does). The properties of the region used as a reference with respect to height are determined by a more global examination of the surface of substrate 60.

An error in the either of the values assumed for reflectivity coefficients R₁ ^(1/2) and R₂ ^(1/2) can introduce an error in the determination of (h₁−h₂). Information about local values of reflectivity coefficients R₁ ^(1/2) and R₂ ^(1/2) can be obtained by an independent measurement. The measured conjugate quadratures in the second embodiment of cited U.S. Provisional Application No. 60/447,254 (ZI-40) and U.S. patent application Ser. No. 10/778,371, filed Feb. 13, 2004 (ZI-40) entitled “Transverse Differential Interferometric Confocal Microscopy” are proportional to components of complex amplitude V₁ (h₁, 0, h₁, 0, χ=π). The complex amplitude V₁ (h₁, 0, h₁, 0, χ=π) of the cited U.S. Provisional Application No. 60/447,254 (ZI-40) and U.S. patent application Ser. No. 10/778,371, filed Feb. 13, 2004 (ZI-40) entitled “Transverse Differential Interferometric Confocal Microscopy” corresponds to a special case of Equation 36, i.e., V ₁(h ₁,0,h ₂,0,χ=π)=[R ₁ ^(1/2) −R ₂ ^(1/2)](1−iβh ₁).  (38) Thus an independent determination the difference in local values of reflectivity coefficients R₁ ^(1/2) and R₂ ^(1/2) can be obtained from a measured value of complex amplitude V₁ (h₁, 0, h₁, 0, χ=π) and used to reduce an error in the determination of (h₁−h₂).

The measurement of differences in reflectivity coefficients can be also used to obtain information about properties of the respective regions of substrate 60 such as information about the complex index of refraction of one the regions relative to the complex index of refraction of a second of the regions.

The selection of values for z₁ in the second embodiment is made to optimize the sensitivity of V₂ (h₁,z₁,h₂,z₂,χ=π) such as expressed by Equations (36) and (37) to quantities for which information is desired in a given end use application. In the selection process, properties of the spherical Bessel functions will also play an important role. One property of j₁ (x) is that j₁ (x) exhibits a maximum for x₁≅2.1 and another property is that j₀ (x)=0 for x≅3.15.

Operation in a dark field mode leads to both reduced systematic and statistical errors in the information represented by the arrays of conjugated quadratures and increased throughput. The information may comprise the transverse derivative of a profile of one or more surfaces of substrate 60 in or on substrate 60; one-dimensional, two-dimensional, and three-dimensional transverse differential images of substrate 60; critical dimensions of features or artifacts on or in substrate 60, and the sizes and locations of sub-wavelength defects in or on substrate 60.

The background components of return measurement beams generated by scattering/reflection of measurement beam components by conjugate spots are the same and therefore do not contribute to the electrical interference signals 72. Accordingly, the background components do not contribute to either the average values of the electrical interference signals 72 or to interference terms in electrical interference signals 72 for both the first and second embodiments.

The reduction of statistical error is also a direct consequence of operation in the dark field mode. The contributions of background fields are removed/eliminated in the second embodiment by the superposition of background fields arranged to have the same amplitudes and phase differences of π and not by the subtraction of intensities. As a result of the dark field, the intensity of beam 24 can be increased significantly without saturation of detector 70 and a corresponding reduction in statistical error is achieved.

The increase in throughput is a direct consequence of operating in a dark field mode. The time required to achieve a certain precision in the measured array of conjugated quadratures is reduced by an increase of the intensity of beam 24 that is permitted by operating in the dark field mode. As a result of the dark field, the intensity of beam 24 can be increased significantly without saturation of detector 70.

Also when operating in a dark field mode, a measured conjugated quadratures of fields corresponding to a pair of spots comprising a sub-wavelength artifact in a locally isotropic section of substrate 60 represents information about the sub-wavelength artifact relative to a reference sub-wavelength artifact. The reference sub-wavelength artifact has properties of the locally isotropic section and dimensions similar to those of the artifact. Accordingly, properties measured include information about critical dimensions and location of the sub-wavelength artifact in or on substrate 60.

Also when operating in a dark field mode, a measured conjugated quadratures of fields corresponding to a pair of spots comprising a sub-wavelength defect in a locally isotropic section of substrate 60 represents information about the sub-wavelength defect relative to a reference sub-wavelength defect. The reference sub-wavelength defect has properties of the locally isotropic section and dimensions similar to those of the defect. Accordingly, properties measured include information about dimensions and location of the sub-wavelength defect in or on substrate 60.

The accuracy of the interferometric compensation of background fields is high in the first and second embodiments for several reasons. The high accuracy of interferometric compensation is not dependent on the properties of pinholes in pinhole array beam-splitter 112, e.g., the diameter of a pinhole could change by a factor of 2 for example and/or the shape of a pinhole could change from a round aperture to a square aperture and the level of interferometric compensation for associated background fields would not change. The amplitudes and phases of background fields associated with a first spot of a pair of spots are the same as the amplitudes and phases of background fields associated with a second spot of a pair of spots independent of properties of pinholes in pinhole array beam-splitter 112.

Compensation for the effects of a mismatch of indices of refraction at the interface of substrate 60 with an external a medium, e.g., air, may be compensated in the first and second embodiments by the addition of a thin low index of refraction layer between pinhole array beam-splitter 112 and lens 150 such as described in U.S. Provisional Patent Application No. 60/444,707(ZI-44) entitled “Compensation for Effects of Mismatch in Indices of Refraction at a Substrate-Medium Interface in Confocal and Interferometric Confocal Microscopy” and U.S. patent application Ser. No. 10/771,785, filed Feb. 4, 2004 (ZI-44) and also entitled “Compensation for Effects of Mismatch in Indices of Refraction at a Substrate-Medium Interface in Confocal and Interferometric Confocal Microscopy” both of which are by Henry A. Hill. The contents of the U.S. Provisional Patent Application and the U.S. Patent Application are incorporated herein in their entirety by reference. With the compensation for effects of the mismatch in refractive indices, an interior portion of a substrate is imaged with a lateral resolution down to of the order of 100 nm and a longitudinal resolution down to of the order of 200 nm. The images of the interior portion are obtained with a working distance of the order of a mm and for depths within the substrate of the order of at least 3 microns.

The throughputs of the first and second embodiments can be further increased by the use of a pinhole array beam-splitter that is coupled to input beam 24 by a guided wave structure such as described in commonly owned U.S. Provisional Patent Application No. 60/445,739 (ZI-39) entitled “Multiple-Source Arrays Fed By Guided Wave Structures And Resonant Structures For Confocal And Near-Field Confocal Microscopy” and U.S. patent application Ser. No. 10/774,250, filed Feb. 6, 2004 (ZI-39) and also entitled “Multiple-Source Arrays Fed By Guided Wave Structures And Resonant Structures For Confocal And Near-Field Confocal Microscopy” both of which are by Henry A. Hill. The contents of the cited U.S. Provisional Patent Application and the U.S. Patent Application are incorporated herein in their entirety by reference.

The first and second embodiments may also be configured for quad-homodyne detection such as described herein and in cited U.S. Provisional Patent Application No. 60/442,858 (ZI-47) and cited U.S. patent application Ser. No. 10/765,369, filed Jan. 27, 2004 (ZI-47) entitled “Apparatus and Method for Joint Measurements of Conjugated Quadratures of Fields of Reflected/Scattered Beams by an Object in Interferometry” (ZI-47).

A third embodiment is shown schematically in FIG. 3. The third embodiment can be configured to be functionally equivalent to the first and second embodiments. The primary difference between the first and second embodiments and the third embodiment is replacement of the pinhole array beam-splitter 112 with traditional confocal pinhole arrays 112A, 112B, and 112C.

Referring to FIG. 3, input beam 24 is incident on polarizing beam-splitter 330 and a first portion thereof is transmitted as a measurement beam of interferometer 410 and a second portion thereof is reflected as a reference beam of interferometer 410 after reflection by mirrors 332 and 334. The measurement beam and the reference beam are incident on pinhole arrays 112A and 112B, respectively. Pinhole arrays 112A and 112B are each conjugates of pinhole array 112C.

A portion of the reference beam incident on pinhole array 112B is transmitted by beam-splitter 340 and a portion thereof focused by lens 360 to an array of spots on pinhole array 112C.

A portion of the measurement beam incident on pinhole array 112A is transmitted by beam-splitter 340 and first and second portions thereof are focused to arrays of spots on substrate 60. The first portion is focused to a first array of spots after a reflection and transmission by polarizing beam-splitter 342, a double pass through quarter-wave plate 346, reflected by concave mirror 350, and focused by lens 354. The second portion is focused to a second array of spots after a transmission and reflection by polarizing beam-splitter 342, a double pass through quarter-wave plate 348, reflected by convex mirror 352, and focused by lens 354. The description the first and second arrays of spots is the same as description of the corresponding arrays of spots of the first and second embodiments. The relative transverse and longitudinal shifting the first and second arrays of spots is controlled by rotations and longitudinal displacements of concave mirror 350 relative to convex mirror 352. The orientation of optical system 310 which in practice is at an angle of 45 degrees to the plane of FIG. 3 is however shown as oriented in the plane in FIG. 3 in order to simplify the description without limiting the scope or spirit of the present invention.

Portions of the measurement beams that form the first and second arrays of spots are reflected/scattered by substrate 60 as first and second arrays of return measurement beam components, respectively. The first array of return measurement beam components retraces the path of its progenitor array of measurement beam components through imaging system 310 and a portion thereof is focused to an array of spots at pinhole array 112C after reflection by beam-splitter 340. The second array of return measurement beam components retraces the path of its progenitor array of measurement beam components through imaging system 310 and a portion thereof is focused to an array of spots at pinhole array 112C after reflection by beam-splitter 340.

The description of the two arrays of spots at pinhole array 112C is the same as portion of the description given for the corresponding arrays of spots in the first and second embodiments at pinhole array beam-splitter 112 except that the displacements of the spots are not reduced by the factor 1/n. The index of refraction of the medium contiguous to pinhole array 112C is assumed to be 1 although it could be otherwise without departing from the scope and the spirit of the present invention.

Portions of the superimposed array of spots and of the reference beam are transmitted by pinhole array 112C and detected by detector 70 after transmission by analyzer 362 to generate electrical interference signal 72. Analyzer 362 mixes the polarization states of the transmitted portions of the superimposed array of spots and of the reference beam.

The description of input beam 24 is the same as corresponding portions of the description given for input beam 24 of FIG. 1 a with beam-conditioner 22 configured as a two-frequency generator and phase-shifter shown in FIG. 1 b and beam 20 comprising a single frequency component. Input beam 24 comprises two components that have different frequencies and each component has two components of different states of plane polarization. The relative phases of the components of input beam 24 are shifted between different values according to control signal 74 generated by electronic processor and controller 80 as described in the discussion of beam-conditioner 22.

The remaining portion of the description of the third embodiment is the same as corresponding portions given for the first and second embodiments.

In some embodiments, pinhole array beam-splitter 112 may be scanned in a direction opposite to the direction of scan of substrate 60 and with a speed such that the conjugate images of the pinholes of pinhole array beam-splitter 12 stay superimposed with spots on or in substrate 60 that are being imaged. This scanning mode of operation is analogous to the relative motions of reticle stage and a wafer stage of a lithography tool operating in a scanning mode. The issue of traditional critical alignment of conjugate confocal pinholes in a confocal microscopy system is nonexsistent, i.e. the registration of the pinholes generating the array of reference beams and the pinholes generating the array of measurement beams is automatic.

A fourth embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 comprising an interferometric far-field confocal microscope such as described in cited U.S. Pat. No. 5,760,901. In the fourth embodiment, beam-conditioner 22 is configured as the two frequency generator and phase-shifter shown in FIG. 1 b. Embodiments in cited U.S. Pat. No. 5,760,901 are configured to operate in either the reflection or transmission mode. The fourth embodiment has reduced effects of background because of background reduction features of cited U.S. Pat. No. 5,760,901.

A fifth embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 comprising an interferometric far-field confocal microscope such as described in cited U.S. Pat. No. 5,760,901 wherein the phase masks are removed. In the fifth embodiment, beam-conditioner 22 is configured as the two frequency generator and phase-shifter shown in FIG. 1 b. Embodiments in cited U.S. Pat. No. 5,760,901 are configured to operate in either the reflection or transmission mode. The fifth embodiment with the phase masks of embodiments of cited U.S. Pat. No. 5,760,901 removed represent applications of confocal techniques in a basic form.

A sixth embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 comprising an interferometric far-field confocal microscope such as described in cited U.S. Pat. No. 6,480,285 B1. In the sixth embodiment, beam-conditioner 22 is configured as the two-frequency generator and phase-shifter shown in FIG. 1 b. Embodiments in cited U.S. Pat. No. 6,480,285 B1 are configured to operate in either the reflection or transmission mode. The sixth embodiment has reduced effects of background because of background reduction features of cited U.S. Pat. No. 6,480,285 B1.

A seventh embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 comprising an interferometric far-field confocal microscope such as described in cited U.S. Pat. No. 6,480,285 B1 wherein the phase masks are removed. In the fifth embodiment, beam-conditioner 22 is configured as the two-frequency generator and phase-shifter shown in FIG. 1 b. Embodiments in cited U.S. Pat. No. 6,480,285 B1 are configured to operate in either the reflection or transmission mode. The seventh embodiment with the phase masks of embodiments of cited U.S. Pat. No. 6,480,285 B1 removed represent applications of confocal techniques in a basic form.

An eighth embodiment comprises the interferometer system of FIGS. 1 a–1 c with interferometer 10 comprising an interferometric near-field confocal microscope such as described in cited U.S. Pat. No. 6,445,453 (ZI-14). In the eighth embodiment, beam-conditioner 22 is configured as the two-frequency generator and phase-shifter shown in FIG. 1 b. Embodiments in cited U.S. Pat. No. 6,445,453 are configured to operate in either the reflection or transmission mode. The eighth embodiment of cited U.S. Pat. No. 6,445,453 in particular is configured to operate in the transmission mode with the measurement beam separated from the reference beam and incident on the substrate being imaged by a non-confocal imaging system, i.e., the measurement beam at the substrate is not an image of an array of pinholes but an extended spot. Accordingly, the corresponding embodiments of the eighth embodiment represent an application of bi-homodyne detection method in a non-confocal configuration for the measurement beam.

Other embodiments may use the quad-homodyne detection method instead of the bi-homodyne detection method as variants of the embodiments. For the embodiments that are based on the apparatus shown in FIGS. 1 a–1 c, the corresponding variants of the embodiments that use the quad-homodyne detection method use variants of the apparatus shown in FIGS. 1 a–1 c. In the variants of the apparatus such as used in the first embodiment, microscope 220 modified to include a dispersive element such as a direct vision prism and/or a dichroic beam-splitter. When configured with a dichroic beam-splitter, a second detector is further added to the system. Descriptions of the variants of the apparatus are the same as corresponding portions of descriptions given for corresponding systems in cited U.S. Provisional Application No. 60/442,982 (ZI-45) and U.S. patent application Ser. No. 10/765,229, filed Jan. 27, 2004 (ZI-45) entitled “Interferometric Confocal Microscopy Incorporating Pinhole Array Beam-Splitter”.

Variants of embodiments may be configured to use the double-homodyne detection method for generation of non-joint measurements of conjugated quadratures. Input beam 24 of the variants of the embodiments comprise four frequency components and with the design of the dispersion of a direct vision prism and/or a dichroic beam-splitter such as described with respect to embodiments that are configured to use the quad-homodyne detection method and the selection of the four frequencies such that each of the four frequency components of beam 32 are directed to different pixels of detector 70. Four arrays of electrical interference signal values are obtained simultaneously and processed for amplitudes of conjugated quadratures using the procedure described herein for the single-homodyne detection method. 

1. An interferometric confocal microscope for measuring an object, said microscope comprising: a source-side pinhole array; a detector-side pinhole array; and an interferometer that images the array of pinholes of the source-side pinhole array onto a first array of spots located in front of an object plane located near where the object is positioned and onto a second array of spots behind the object plane, wherein the first and second arrays of spots are displaced from each other in both a direction normal to the object plane and a direction parallel to the object plane, said interferometer also imaging the first arrays of spots onto a first image plane that is behind the detector-side pinhole array and imaging the second array of spots onto a second image plane that is in front of the detector-side pinhole array wherein each spot of the imaged first array of spots is aligned with a corresponding different spot of the imaged second array of spots and a corresponding different pinhole of the detector-side pinhole array.
 2. A differential interferometric confocal microscope for measuring an object, said microscope comprising: a source-side pinhole array; a detector-side pinhole array; and an interferometer that images each pinhole of the source-side pinhole array onto a corresponding different pair of two locations, one of which lies in a first object plane and the other of which lies in a second object plane that is parallel to and displaced from the first object plane, thereby generating a first image of the source-side pinhole array in the first object plane and a second image of the source-side pinhole array in the second object plane, said interferometer also projecting a first array of return measurement beams from the first image and a second array of return measurement beams from the second image toward the detector-side pinhole array to produce a first array of converging beams and a second array of converging beams, wherein the detector-side pinhole array generates an array of fields each of which is a difference of fields of the first and second arrays of converging beams.
 3. A differential interferometric confocal microscope for measuring an object and which has, in the vicinity of where the object being measured is to be located, a first object plane and a second object plane that is displaced from and parallel to the first object plane, said microscope comprising: a source-side pinhole array; a detector-side pinhole array; and an interferometer that receives a beam from a selected pinhole of the source-side pinhole array and converges a first part of that received beam onto a corresponding first location in the first object plane and a second part of that received beam onto a corresponding second location in the second object plane, said interferometer also arranged to receive a first return beam from the first location and a second return beam from the second-location and converge at least a part of each of the first and second return beams onto a corresponding pinhole of the detector-side pinhole array to produce a difference of fields of the first and second return beams converging on that corresponding pinhole, wherein said selected pinhole is any pinhole of the source-side pinhole array.
 4. An interferometric confocal microscope for measuring an object, said microscope comprising: a source-side pinhole array for producing an array of input beams; a detector-side pinhole array; and an interferometer including: a first optical element providing a first reflecting surface; a second optical element providing a second reflecting surface; and a beam splitter positioned between the first and second optical elements, wherein the beam splitter produces from the array of input beams a first array of measurement beams and a second array of measurement beams, wherein the first reflecting surface participates in focusing the first array of measurement beams onto a first array of locations on a first object plane in object space and the second reflecting surface participates in focusing the second array of measurement beams onto a second array of locations on a second object plane in object space, said first and second object planes being parallel to and displaced from each other, wherein the first array of measurement beams generates a first array of return beams from the object and the second array of measurement beams generates a second array of return beams from the object, wherein the first reflecting element participates in producing from the first array of return beams a first array of converging beams that converge to a first array of spots on a first image plane and the second reflecting element participates in producing from the second array of return beams a second array of converging beams that converge onto a second array of spots on a second image plane, said first and second image planes being adjacent to and on opposite sides of the detector-side pinhole array, and wherein the detector-side pinhole array combines the first and second arrays of converging beams to form an array of output beams.
 5. The interferometric confocal microscope of claim 4 wherein a single pinhole array serves as both the source-side pinhole array and the detector-side pinhole array.
 6. The interferometric confocal microscope of claim 5, wherein the first optical element is located between said single pinhole array and the beam splitter and wherein the second optical element is located between a location at which the object is positioned during use and the beam splitter, wherein the first reflecting surface has a center of curvature for which there is a corresponding conjugate as viewed through the beam splitter, and wherein the second reflecting surface has a center of curvature that is displaced relative to the corresponding conjugate of the center of curvature of the first reflecting surface.
 7. The interferometric confocal microscope of claim 6, wherein the conjugate of the center of curvature of the first reflecting surface and the center of curvature of the second reflecting surface are displaced from each other in a first direction that is normal to a plane defined by the beam splitter and in a second direction that is parallel to the plane defined by the beam splitter.
 8. The interferometric confocal microscope of claim 7, wherein the first reflecting surface participates in focusing the first array of measurement beams via the beam splitter onto the first array of locations and the second reflecting surface participates in focusing the second array of measurement beams via the beam splitter onto the second array of locations.
 9. The interferometric confocal microscope of claim 8, wherein the first reflecting element participates in combination with the beam splitter in producing the first array of converging beams and the second reflecting element participates in combination with the beam splitter in producing the second array of converging beams.
 10. The interferometric confocal microscope of claim 9 wherein the first reflecting surface is substantially concentric with a point on the object.
 11. The interferometric confocal microscope of claim 9, wherein the second optical element provides a refracting surface positioned between the object and the beam splitter to receive light rays from the object.
 12. The interferometric confocal microscope of claim 11, wherein the first reflecting surface substantially conforms to a sphere having a first radius and the refracting surface conforms to a sphere having a second radius, wherein the first radius is greater than the second radius.
 13. The interferometric confocal microscope of claim 9, wherein the first optical element provides a refracting surface positioned between the beam splitter and said single pinhole array.
 14. The interferometric confocal microscope of claim 9 wherein the second reflecting surface is substantially concentric with an image point on said single pinhole array.
 15. The interferometric confocal microscope of claim 13, wherein the second reflecting surface substantially conforms to a sphere having a first radius and the refracting surface conforms to a sphere having a second radius, wherein the first radius is greater than the second radius.
 16. The interferometric confocal microscope of claim 9, wherein said single pinhole array is a two-dimensional array.
 17. The interferometric confocal microscope of claim 16, wherein the two-dimensional array is of equally-spaced holes.
 18. The interferometric confocal microscope of claim 17, wherein the equally-spaced holes are circular apertures.
 19. The interferometric confocal microscope of claim 9, wherein the first and second object planes are separated from each other on the order of the longitudinal resolution of the confocal interferometric microscope.
 20. The interferometric confocal microscope of claim 7, wherein the conjugate of the center of curvature of the first reflecting surface and the center of curvature of the second reflecting surface are displaced by an amount that causes the detector-side pinhole array to generate an array of fields each of which is a difference of fields of corresponding beams of the first and second arrays of return measurement beams.
 21. The interferometric confocal microscope of claim 7, wherein the conjugate of the center of curvature of the first reflecting surface and the center of curvature of the second reflecting surface are displaced by an amount that causes the interferometric confocal microscope to operate in a differential mode. 